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Presentacion 



Desde que Foynman en 1948 introdujera la integral de camino en Teoria 
Cuantica de Campos (TCC), una gran parte de los fisicos teoricos de 
particulas han trabajado usando este formalismo. 

Las predicciones obtenidas estuvieron practicamente limitadas al regi- 
men perturbativo (pequenas fluctuacioncs cn torno al campo libre) en los 
primeros anos, puesto que en este limite las integrales son Gaussianas y por 
tanto se pueden resolver por metodos analiticos. 

No perturbativamente, una posible regularizacion se obtiene reemplazando 
el espacio-tiempo continue por una red discreta (Wilson, 1974). Por supuesto 
la prcgunta a responder es si las observaciones fisicas que tenemos en el con- 
tinue se corresponden de alguna mancra con la formulacion de la teoria en 
la red. Si se puede probar que el limite continue de una teoria de campos 
en la red existe, uno podria usar la definicion de la teoria en la red como 
definicion de la teoria en el continue. 

El trabajo pionero de Wilson desencadeno interes por atacar la solucion 
de las integrales funcionales numericamente. El impresionante auge de la 
informatica cn la ultima decada, con el consiguiente aumento de los recursos 
accesibles de potencia de calculo, ha favorecido sin duda el avance de este 
tipo de estudios, gracias a los cuales se ha obtenido conocimiento sobre los 
aspectos no pertubativos de TCC. 

A lo largo de esta memoria se van a estudiar las propiedades de distintos 
modelos relevantes para la Fisica de Particulas mediante su formulacion 
en una red espacio-temporal. Los analisis estaran basados en resultados 
que se obtienen con la ayuda de ordenadores. Para ello se van a usar 
ideas provenientes de la Mecanica Estadistica Clasica. Asi utilizaremos el 
concepto de transicion de fase y de parametro de orden, aplicandolos en 
diferentes contextos. 

Cualquiera que este interesado en el estudio numerico de fenomenos 
crfticos dcbe considcrar como fundamental el hecho de quo las transicioncs 
de fase ocurren solo en el limite termodinamico. Cuando se hacen si- 
mulaciones numericas se esta forzado a trabajar con un mimero finito de 



grades de libertad, por tanto no hay transiciones de fase. Sin embargo 
los efectos que se observan en redes de tamaiio finito son precursores del 
verdadero comportamiento Hmite. La teoria de Finite Size Scaling dicta la 
forma concreta en que esto ocurre y por lo tanto, mediante el estudio de la 
evolucion de los observables de interes con el tamaiio de la red se pueden 
obtener predicciones sobre el comportamiento en el Ifmite termodinamico. 
Las tecnicas de tamaiio finito son por lo tanto una herramienta fundamental 
en los analisis que se realizan a lo largo de esta memoria. 

La estructura de la memoria es la siguiente. Se ha querido, por mo- 
tivos de completitud, hacer una breve introduccion al estudio de fenomenos 
crfticos en la red en un capitulo preliminar. Tras este preliminar, en primer 
lugar se describe el modelo 0(4) Anti-Ferromagnetico. La motivacion ini- 
cial de este trabajo fue la busqueda de puntos fijos no triviales en cuatro 
dimensiones en el sector Anti-Ferromagnetico del Modelo Estandar. 

El segundo capftulo corresponde al estudio de una transicion que podriamos 
decir clasica en el marco do las simulaciones en la red de TCC, la transicion 
de fase en el modelo SU(2)-Higgs a temperatura cero. 

En el tercer capftulo se ha abordado un problema que durante los 
liltimos afios esta sometido a una vigorosa discusion: la influencia de las 
condiciones de contorno en el orden de la transicion de fase de la teorfa 
U(l) compacta pura gauge. 

Debido a la conexion existente entre Mecanica Estadfstica y TCC es - 
natural realizar incursiones en problemas puramentc dc Mecanica Estadfstica, 
dentro de la interaccion con otros grupos de invcstigacion. En mi caso esta 
incursion se hizo estudiando un problema podrfamos decir que de moda 
como es la formulacion de modelos adaptados a la descripcion de proce- 
sos que envuelven un flujo de informacion. Este trabajo se describe en el 
capftulo final de esta memoria. 

La estructura de cada capftulo consta de una introduccion para motivar 
y presentar el trabajo realizado, seguida de una exposicion del metodo 
utilizado y de los resultados obtenidos. 

Finalmente se comentan las conclusiones obtenidas y se detalla la lista 
de publicaciones. 
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Chapter 1 

Preliminares 
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1.1 Regularizacion no perturbativa en Teoria 
Cuantica de Campos 

La Teoria Cuantica de Campos es el marco mas adecuado para describir 
las interacciones fuerte y electrodebil. En el Modelo Estandar (SM) ^ la 
teoria unificada que describe estas interacciones esta basada en el grupo 
gauge SU{3) (g> SU(2) ® U{1)y. 

Uno de los posibles formalismos que permiten estudiar las propiedades 
de una TCC es el de la integral de camino. A esta formulacion se llega a 
traves de la generalizacion del concepto de integral de camino introducido 
per Feynmann para la Mecanica Cuantica, y su extension a tiempos imagi- 
narios (integral de camino Euclidea). En este contexto las funciones de 
Green Eucli'deas o funciones de Schwinger a partir de las cuales se construye 
la TCC [||, se pueden escribir como los momentos de una determinada 
medida de probabilidad: 

5„(a;i,...,a;„)- 1 [d$](f>(xi) • • • <I)(x„)e-^[*l , (1.1) 

donde 5'[$] es la accion Euclidea 

= j dtC{^,df,<^) (1.2) 

La formulacion de la integral de camino no solo proporciona una vision 
ffsica llamativa de la evolucion de un sistcma cuantico como una suma 
sobre caminos clasicos, sino que abre el camino a calculos tanto analfticos 
como numericos. Ademas ayuda a comprender la relacion entre TCC y 
fenomenos criticos en Mecanica Estadistica (ME) clasica. Para ver esta 
relacion consideremos la funcion de correlacion a n puntos de una TCC 
escalar en el espacio Euclideo: 

(0|r(<i>(xi) • • • $(x„))|0) = I y [d^Mx^) ■ ■ ■ $(x„)e-^[*l , (1.3) 

donde hemos incluido la normalizacion Z de manera que el valor medio de 
la identidad sea uno: 

f [d$]e-^[*l . (1.4) 



Estas ecuaciones se pueden interpretar como la funcion de correlacion y 
la funcion de particion respectivamente de un sistema de ME clasica. La 
TCC en ds dimensiones espaciales se puede hacer equivalente a un sistema 
de ME clasica en equilibrio en d = dg + 1 dimensiones Euclideas. La accion 
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EucUdea juega el papel de E/kbT, donde con E se denota la energia clasica 
de la configuracion Q 

Una primera aproximacion a la solucion de este tipo de integrales fun- 
cionales es considerar pequenas fluctuaciones en torno al campo libre (in- 
tegral gausiana), este es el contexto de la Teoria de Perturbaciones. La 
resolucion perturbativa pasa por hacer un desarroUo en serie de la parte 



no gausiana de (1.3). Sin embargo al calcular la contribucion de los distin- 
tos terminos de la serie perturbativa aparecen divergencias asociadas a los 
grandes momentos como consecuencia de trabajar con un numero infinito de 
grados de libertad (la variable x es continua). En los calculos perturbativos 
estas divergencias ultravioletas (UV) han de ser compensadas mediante la 
introduccion de contratcrminos en el lagrangiano de partida, lagrangiano 
desnudo, expresado en funcion de parametros desnudos, acoplamientos y 
masas. 

Para hacer calculos hay que dar sentido a estas integrales divergentes, 
mediante un cut-ofF en momentos, trabajando en dimension menor que 4, 
etc. Este procedimiento se denomina regularizacion. Una vez que se han 
hecho los calculos se suprime la regularizacion, para ello se fijan un cierto 
numero de parametros (acoplamientos y masas) a sus valores renorma- 
lizados y se hace tender el cut-off a infinito o la dimension a 4. En este 
proceso, los infinitos que aparecen al calcular los terminos de la serie per- 
turbativa son absorbidos en los parametros desnudos y las constantes de 
renormalizacion que aparecen en los contraterminos. Una teon'a se dice 
renormalizable cuando bastan un numero finito de contraterminos para 
hacer la teon'a finita a todos los ordenes de teorfa de perturbaciones. En 
dimension 4 las teorfas gauge, como el SM, son renormalizables y por 
tanto no hay ningiin problema en dar sentido perturbativo a la dcfinicion 



(1.1) 



Resumiendo, las integrales funcionales del tipo (1.1) no estan definidas 
apropiadamente, su expresion es puramente formal y para darle sentido hay 
que regularizar la teorfa de manera que no aparezcan infinitos en el UV. 
En teorfa de perturbaciones son posibles varias regularizaciones que nos 
permiten dar sentido a las integrales. 

Sin embargo, si se quiere dar cuenta de fenomenos no accesibles al 
calculo perturbativo, como son por ejemplo la rotura espontanea de la 
simetrfa quiral, o la generacion de masas para los hadrones, necesitamos 
una regularizacion que permita dar sentido al formalismo de la integral 
de camino en el regimen no perturbativo. Se obtiene una regularizacion 
no perturbativa introduciendo una longitud elemental "a" en el espacio- 
tiempo En este esquema las coordenadas del espacio-tiempo estan 



^ Se ha usado el convenio fc;, = 1 a lo largo de esta memoria 
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discretizadas: 

Xfj. = aUf,. Uf, e Z . (1.5) 
Debido a la discretizacion, la transformada do Fourier de una funcion G(x) 

G{k) = a^Yl e'"''^"''G(a;) , (1.6) 

es periodica, es decir invariante bajo 

kf,^kf, + 27r/a . (1.7) 

For lo tanto se puede restringir fc^ a todo intervalo de longitud 27r/a, por 
ejemplo — tt/o < fc^ < ""/oi) Que recibe el nombre de primera zona de 
Brillouin. Es decir, la consccucncia de haber introducido una longitud 
elemental "a" , es que hemos regularizado la teoria en el UV pucsto que 
ahora los momentos estan restringidos a los valores del interior dc una caja 
de arista 2n/a. Esta rcgularizacion es mas potente que las regularizaciones 
perturbativas, puesto que se espera que sirva como definicion de la teorfa 
completa, no solo de su aspecto perturbativo. 

1.2 Formulacion en la red de la TCC 

Para una teorfa oscalar hemos visto que si colocamos las variables dinamicas 
en los nodos dc la red, el sistema es cquivalente a un modelo de ME donde 
las variables son los valores de los campos. 

En general en los nodos de la red se colocan los campos cscalares y los 
fermionicos, las derivadas se sustituyen por difcrcncias finitas y los campos 
vectoriales se colocan en los segmentos que uncn nodos vecinos (links). 
En cuanto a estos liltimos, parecerfa que la eleccion natural es colocar los 
campos vectoriales (potenciales gauge) en las links. Sin embargo hacerlo asi 
rompe la invariancia gauge poniendo en peligro la renormalizabilidad. La 
forma de poner los campos gauge en la red es poner en las links no el mismo 
campo gauge, sino los clcmcntos del grupo gauge. La link que conecta el 
site n con el site de coordcnadas n + /t se denota por [/^(n), 

C/^(n)=e*«^'<M^'' , (1.8) 

donde fi es el vector unitario en la direcion espacio-temporal /x, y 7^ son los 
generadores del grupo gauge que estemos considerando. 

En la red cs convcnicnte trabajar con magnitudes adimensionales de 
manera que los operadores definidos en el continuo se redefinen en la red 
de la siguiente forma: 

Qcont ~* Qlat = d^Qcont j (1-9) 
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donde s es la dimension de Qcont en unidades de energia. 

La formulacion de la integral de camino en TCC y su regularizacion 
en la red, permite la evaluacion numerica de las integrales funcionales. 
En particular estaremos interesados en calcular el promedio de operadores 
invariantes gauge: 

(0($,A)> = [rf$][dA]0($,A)e-^[*''^l , (1.10) 

en una red finita cuadridimensional de lado L mediante integracion numerica. 



Para obtener informacion sobre (1.10) numericamente definimos 



0= lim / dtOmt),A{t))/Q , (1.11) 

Jo 



donde t es un tiempo de ordenador, definido mediante la dinamica que 
nosotros hayamos elegido. 

En estas condiciones, se evalua el operador O y se suma sobre las dife- 
rentes configuraciones de los campos [$(t), A(r)] generadas en el curso de 
la dinamica, que en general es discreta, es decir, se g eneran configuraciones 
numeradas 1,2,3 etc., con lo cual la integral (|l . 1 1| ) sera en realidad una 
suma del operador O evaluado sobre las diferentes configuraciones. 

La dinamica se suele elegir definiendo una matriz de probabilidad de 
transicion de una configuracion a otra, que sea ergodica (para que el sistema 
no quede atrapado en un subconjunto del espacio fase), y que deje invariante 
la distribucion de probabilidad de Boltzmann. Una tal dinamica se llama 
proceso de Markov y la secuencia de configuraciones generadas se llama 
cadena de Markov. Es posible demostrar que una dinamica que cumple 
estas condiciones converge a la distribucion de equilibrio de Boltzmann 



definida en ( 1.10 ), y por tanto se tiene que: 

d^{0) (1.12) 

En la practica hay dos problemas a tener en cuenta. El primero es que no se 
puede generar un mimero infinito de configuraciones. Cuando se considera 
un mimero finito de configuraciones N, la expresion ( 1.12| ) es cierta con un 



error que es del orden de 1/\N (teorema del Ifmite central). El segundo 
problema es que la dinamica que usamos genera configuraciones que estan 
correlacionadas. El tiempo que tarda el sistema en perder la memoria de la 
configuracion de la que viene se denomina tiempo de autocorrelacion y se 
denota por r. Hablando cualitativamente, si se generan N configuraciones 
el mimero de configuraciones independientes es N /t. Para una exposicion 
detallada sobre los errores intrfnsecos a una simulacion de Monte Carlo 
vease [H. 
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Del mismo modo que se hace en teon'a de perturbaciones, una vez que 
hemos dado sentido a las integrales hay que eliminar la regularizacion. El 
proceso de eliminar la regularizacion impuesta per la formulacion en la 
red se denomina tomar el li'mite continue. La forma precisa en que hay que 
tomar este h'mite la dicta el Grupo de Renormalizacion (RG) . En particular 
cuando se esta suficientemente proximo al h'mite continue las magnitudes 
ffsicas no deben depender de los acoplamientos desnudos ni de "a" . Es decir, 
la renormalizabilidad de la teorfa implica la existencia de trayectorias de 
fisica constante en el espacio de acoplamientos desnudos. 

En este contexto el proceso a seguir sen'a el siguiente. Se evalua en la red 
el operardor Q{g) correspondiente al observable Q para distintos valores de 
los acoplamientos g. Despues se restauran las dimensiones multiplicando 
Q{g) por la potencia adecuada de "a", obtenemos asf Q{g,a). La forma 
en que depende g de "a" la sabemos perturbativamente, de manera que al 
menos en la region perturbativa sabemos por donde pasan las trayectorias 
de fisica constante. Usando pues la prediccion perturbativa podemos cono- 
cer Q{g{a),a). Si el observable Q no depende de "a" podremos decir que 
estamos estudiando correctamente la fi'sica del continue mediante nuestra 
formulacion en la red. Esto es le que se cenece come scaling asintdtico. Una 
condicion mas suave que la anterior es lo que se conoce come scaling. Aquf 
se exige que cecientes adimensionales de cantidades con dimensiones (por 
ejemplo mQ++ /m2++) permanezcan constantes al variar los parametres y 
acercarnos al h'mite continue. Por supueste scaling asintotico implica scal- 
ing pero no a la inversa. 

Supongamos que queremos calcular el espectro de masas de la teorfa, 
calculando para ello el decaimiento en la red de las funciones de correlacion 
aprepiadas para largos tiempos Euch'dees. La masa mas pequefia de la 
teorfa esta determinada por la longitud de correlacion ^ mas larga, que se 
define como la escala de longitud sebre la cual la partfcula se puede propagar 
con una amplitud significativa. La renormalizabilidad de la teorfa exige que 
la masa fi'sica sea finita, por lo tante la masa en unidades de la red, rh = am, 
tiene que ser cere. Esto implica que la longitud de correlacion medida en 
unidades de la red ^ debe diverger. 

Resumiendo, el Ifmite continue de una teorfa de campes en la red puede 
tener lugar sole en un punte crftico del sistema de ME al cual es equivalente 
la TCC, y que esta definide mediante la funcion de particion (L4). Esta 
condicion es por otra parte logica, pueste que sole cuando f — + cx3 el sistema 
pierde memoria de la red subyacente. 
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1.3 Transiciones de fase en el reticulo 



Los mecanismos de transicion de fase aparecen frecuentemente en fi'sica de 
particulas para dar cuenta de las propiedades observadas en la naturaleza. 
Asi por ejemplo, la generacion de masa para los bosones vectoriales se 
entiende a traves del mecanismo de Higgs, que se origina en una transicion 
de fase producida por la ruptura espontanea de la simetri'a del campo de 
Higgs. 

La ausencia de una transicion de fase en SU(2) y SU(3) entre la region 
de acoplamiento fuerte y acoplamiento debil, sugerida por los trabajos pi- 
oneros de Creutz Q, se toma como criterio para pensar que la propiedad 
de confinamiento que puede ser demostrada en el regimen de acoplamiento 
fuerte ||], se mantiene en el regimen de acoplamiento debil, y por tanto 
dan fundamento a la creencia de que la QCD es una teorfa razonable para 
estudiar la interaccion fuerte. 

Como se ha apuntado anteriormente la construccion de una teorfa en el 
contfnuo, implica la existencia de un punto crftico en el sentido de la teori'a 
de transiciones de fase en ME. 

Por lo tanto es fundamental el estudio del diagrama de fases de las TCC 
en el retfculo, tanto desde el punto de vista de la caracterizacion de las fases 
en que se manifiesta la teoria, como para encontrar puntos criticos en los 
cuales la teoria es susceptible de definir una TCC en el continuo. 

Las propiedades de un sistema fisico dependen de parametros como la 
temperatura, la presion, el campo magnetico aplicado, etc... Sin embargo 
un cambio en las condiciones a que esta sometido puede hacer que las 
propiedades del sistema vari'en bruscamente. Este cambio brusco se denomi- 
na transicion de fase. Un cambio brusco, en el lenguaje de las matematicas 
es una discontinuidad, la cual se manifestara en las medidas de las funciones 
termodinamicas. 

Vamos a centrarnos en las transiciones inducidas por la competicion 
orden-desorden. Para tener un soporte concreto la discusion se hara para el 
modelo de Ising. A pesar de su simplicidad las conclusiones que se obtengan 
seran aplicables con gran generalidad. La accion del modelo de Ising en 
presencia de un campo magnetico externo es: 

5 = -/3^a,a, -/i^a. , (1.13) 

donde (3 — J/T y h — H/T, siendo J la intensidad del acoplamiento 
entre espines y H el campo magnetico externo aplicado. El espi'n a toma 
los valores +1 6 -1. La notacion indica que la suma se extiende a las 
parejas de primeros vecinos. 
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La funcion de particion del sistema esta dada por: 

Zat = ^ e-^(^°"f^ (1.14) 

conf 

donde el sumatorio en configuraciones significa: 

E - E E ••• E ' (1-15) 

conf cri=±lcr2=±l criv=±l 

sicndo cl niimcro total de espines. Por tanto el niimero de terminos 
en la funcion de particion es 2^. 

La densidad de energia libre esta dada por: 

F{/3,h) = -^\ogZN , (1.16) 

cuyas dcrivadas con respecto a /3 y/6 a /i nos daran las magnitudes ter- 
modinamicas. 

Definimos el valor medio del espm o-j como: 

(^i) = ^E^'^"^^™"'^ • (1-1'^) 

conf 

Debido a la invarianza traslacional de la accion, y trabajando en una red 
con condiciones dc contorno pcriodicas, este promedio no depende del espm 
i particular para el que se calcula. 

La magnetizacion del sistema estara dada por el valor medio de la suma 
de los espines individuales: 

M{T,h) = N-\Y,ai) , (1.18) 

i 

Denotando por w{E) cl numcro do configuraciones con cncrgia interna 
E, podemos reagrupar la funcion dc particion, y expresarla como una suma 
sobre energias en lugar de como una suma sobre configuraciones: 

Z = J2 w{E)e-^ = E e^[^-™(^)l , (1.19) 

E E 

donde hemos usado la dcfinicion dc cntropfa, Vl = log u>(i<^). La probabili- 
dad de aparicion de una configuracion con cncrgfa interna E es proporcional 
a su factor de Boltzman multiplicado por la dcgcncracion w{E). Si el estado 
fundamental tiene una degeneracion finita, la configuracion corrospondicntc 
a T = sera aquella que corresponda al estado de minima energia (todos los 
espines alineados) pues es infinitamente mas probable que cualquier otra. 
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El modelo de Ising presenta una Imea de transiciones de fase a campo 
magnetico cero: {h = 0,0 < T < Tc}. Que las transiciones de fase tienen 
que ocurrir a campo magnetico nulo es una consecuencia de la simetrfa de 
la accion ||] : 

a^^-(Ti, h^-h. (1.20) 

Por encima de Tc las fluctuaciones termicas desordenan por completo al 
sistema y ya no hay transiciones de fase. 

El parametro conveniente para estudiar la aparicion de la transicion de 
fase es la magnetizacion, en particular la magnetizacion espontanea definida 
como: 

Mo(T) = lim M(T, h) . (1.21) 

Por encima de Tc la magnetizacion espontanea se anula, mientras que 
es distinta de cero para T <Tc- 

Las transiciones de fase se clasifican de acuerdo con el comportamiento 
cualitativo de Mq en el punto de la transicion. En el modelo de Ising Mq 
es continua en T = Tc, y acorde a ello la transicion se dice que es continua 
en el punto {h^O,T ^Tc}. 

Sin embargo, en el resto de puntos de transicion de fase de la Ifnea, es 
decir {h — 0,0 <T < Tc\, Mq se comporta de forma discontinua: 

lim M(T,h) = ~Mo(T) 

h^O- 

\im M(T,h) = +Mo(T) . (1.22) 

Correspondientemente se dice que las transiciones de fase son de primer 
or den o discontinuas. 

Por ser la magnetizacion espontanea la magnitud que indica en que 
fase se encuentra el sistema, recibe el nombre de parametro de orden de la 
transicion. 

Hablando con generalidad, se dice que una transicion de fase es de 
primer orden cuando hay una discontinuidad en la primera derivada del 
potencial termodinamico, en el caso que estamos considerando, de la en- 
ergi'a libre: 

Energfa interna: U = dF/dB 

(1.23) 

Magnetizacion: M — —{dF/dh) 

Si las primeras derivadas son continuas, pero aparecen discontinuidades 
en las derivadas segundas (al menos una de ellas es discontinua) la transicion 
se dice de segundo orden. En general se habla de transiciones de fase 
continuas cuando la discontinuidad aparece en una derivada de orden igual 
o superior a dos, siendo todas las primeras derivadas continuas. 
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En el caso Ising en el punto {h ^ 0,T ~ Tc) la transicion es de se- 
gundo orden por lo que aparecen divergencias en el calor especi'fico y en la 
susceptibilidad magnetica: 

Calor Especi'fico: = dU/d(3 = (d'^F/df3^) 

(1.24) 

Susceptibilidad: x ^ dM/dh = -{d^F/dh'^) 

La diferencia fundamental entre ambos tipos de transiciones es que las tran- 
siciones continuas presentan ordenamientos de largo alcance en el punto de 
la transicion. Las divergencias que aparecen estan ligadas a la divergencia 
de la longitud de correlacion. Sin embargo en las transiciones de primer 
orden el sistema no anticipa la transicion y la longitud de correlacion per- 
manece finita cuando nos aproximamos al punto de la discontinuidad. 

Desde el punto de vista de la TCC las transiciones interesantes son por 
tanto las de segundo orden pues solo en estos puntos se puede obtener una 
teorfa estrictamente renormalizable. Puesto que en esta memoria estaremos 
especialmente interesados en la diferenciacion de transiciones de fase de 
primer y segundo orden, veamos mas detalladamente las caracterfsticas de 
ambas. 



1.3.1 Transiciones de fase continuas 

Para una descripcion cualitativa de los fenomcnos que ocurren en estas 



transiciones de fase tomemos como ejemplo el modelo de Ising ( 1.13 ) en 
ausencia de campo magnetico externo (/i = 0). 

A temperatura elevada se observan espines hacia arriba y espines hacia 
abajo formando islotes cuyo tamario promedio es la longitud de correlacion 
entre espines. La correlacion entre dos espines situados a una distancia 
tiende a cero exponencialmente con la separacion: 

^ e-'--/«(^) , (1.25) 



con una longitud de correlacion caracterfstica de la temperatura, ^(T). En 
estas condiciones la magnetizacion M es nula. 

A medida que disminuye T, el tamafio de los islotes aumenta pues la cor- 
relacion entre espines aumenta al disminuir la agitacion termica. Cuando se 
alcanza la temperatura crftica, T = Tc los islotes tienen todos los tamaiios 
posibles. Es decir, a Tc las fluctuaciones del sistema son de todas las lon- 
gitudes, ya no hay una escala tfpica para el sistema, por eso se dice que 
la fisica es invariante de escala en el punto de la transicion, que recibe el 
nombre de punto cn'tico. 
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Asi mismo por debajo de Tc la funcion de correlacion entre espines 
situados a distancia rij es el parametro de orden al cuadrado: 

^]imjaiaj) = {ai){aj) = . (1.26) 

Se define la funcion de correlacion entre dos espines de la siguiente 
forma: 

r{e,rij) = {aiaj)-{ai){aj) , (1.27) 

donde e = (T — Tc)/Tc recibe el nombre de temperatura reducida. 

Con esta definicion para T > Tc y T < Tc , T{e,rij) decrece exponen- 
cialnicntc con la distancia: T{e,rij) ~ exp(— rjj/^), siendo ^ la longitud 
de correlacion asociada a la funcion de correlacion entre dos espines. La 
definicion precisa viene dada por: 

^ = lim , "'^'^l ^ . (1.28) 
|r„Hoo logr(e,roi) 

Que cxista una magnctizacion cspontanca a campo exterior nulo, es de 
por si un hecho remarcable puesto que uno no esperarfa en estas condiciones 
que el sistema elija una direccion privilegiada hacia la que apuntar. Sin 
embargo la mas pcquciia inlioincgcncidad en los islotes. o un campo residual 
no nulo pueden hacer que una direccion resulte privilegiada frente a las 
demas, o dicho de otra forma el vacfo con M = no es estable. Este 
fenomeno se llama ruptura espontanea de la simetria. 

El hecho de que la transicion se produzca como un fenomeno coopera- 
tivo a gran escala induce a pcnsar que ciertas de sus caracteristicas solo 
dcpcndcn de las propicdadcs del sistema a gran escala. osto os, do sus 
propiedades generales, tales como la dimension d del espacio, la dimension 
del parametro de orden o las simetrfas de los acoplamientos, y no de los de- 
talles de la interaccion. Esta es la idea sobrc la que dcscansa el conccpto de 
Universalidad, que nos llevara a definir magnitudes para describir el com- 
portamiento crftico de los sistemas dependiendo solo de estas caracteristicas 
generales, por ejemplo los exponentes crfticos. 

Los exponentes crfticos caracterizan el comportamiento de los parame- 
tros de orden, susceptibilidades, etc, en el entorno del punto crftico. La 
hipotesis que se hace es que el comportamiento de los observables en el 
entorno del punto crftico se puede describir como una potencia de la tem- 
peratura reducida, e. Esta potencia se llama exponente crftico y cada ob- 
servable ffsico dcpcndicntc dc la temperatura tiene asociado uno. De esta 
manera dado un observable /(e), sc hace la hipotesis de que /(e) es contfnuo 
y positivo para valores pequefios y positivos del parametro e. Se supone 
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ademas que el limite 

A = lim ^^iM , (1.29) 
£^0 loge 

existe y esta bien definido. A A se le llama exponente critico asociado a 

Para sistemas magneticos como el modelo de Ising, los exponentes se 
definen de la siguiente forma: 

a(r) ~ e-" ; 

M(T) - (-e)'3 • 

xiT) ^ 



(1.30) 



Por ultimo se definen tambien exponentes para el comportamiento en 
el punto critico (e=0) dc la funcion de correlacion de dos espines situados 
a distancia r: 

r(0,r)-r-('^-2+.',) ^ (^ 3^) 

En general la funcion de correlacion se puede expresar como: 

ne,r) = ^. (1.32) 

Para T ^ Tc, g{r,S,) ''^ exp(— r/^) mientras que en el entorno del punto 
critico (^ — oo) la funcion de correlacion se comporta como una potencia 
de r. 

Se introduce tambien un exponente critico asociado al comportamiento 
del parametro de orden en el punto critico a campo cxterno distinto de cero: 

M{T,,h)^\h\'/^ . (1.33) 

En la dimension crftica superior estas predicciones estan modificadas 
por la existencia de correcciones logaritmicas . 

Hay dos motivos esencialmente por los cuales estamos interesados en 
estas cantidades a pesar de que posean menos informacion que la funcion 
completa /(e): 

1. Cerca del punto critico (T « Tc) domina el termino con e elevado a la 
potencia mas baja. Esto se confirma experimentalmente en graficas 
log-log puesto que se obtienen rectas en torno al punto critico cuya 
pendiente es el exponente critico. Por otra parte los fenomenos fisicos 
en los que estamos interesados, tales como las existencia de largas 
correlaciones, ocurren en el entorno del punto critico. 
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2. Existen relaciones entre los exponentes que trascienden a cualquier 
sistema particular pues dependen de caracteristicas muy generales, 
en la li'nea del concepto de Universalidad mencionado anteriormente. 
Estas relaciones no han sido probada con la mayor generalidad sino 
solo en el marco de ciertas hipotesis que conciernen a la forma de los 
potenciales termodinamicos. 

En particular se supone que el potencial de Gibbs es una funcion homo- 
genea generalizada {hipotesis de escala) es decir: 

G(A°^e, A°«i/) = AG(e, H) . (1.34) 

Esta hipotesis permite establecer relaciones entre los exponentes cn'ticos: 

7 = 1/(2 - 77) ; 
a + 2/5 + 7 = 2 ; 

1.35 

7-/3(5-1); 
7 = 2/3(5 + a - 2 . 

Si ademas se asume la hipotesis de hiperescala que consiste en suponer 
que las fluctuaciones que contribuyen a la parte singular de la densidad de 
energia libre van como: 

Fsin^ie) ^ ^ ~ (1.36) 

teniendo en cuenta que el calor especi'fico es la segunda derivada de la 
energia libre, podemos escribir su ley de escalado como: 

- e"""^-^ . (1.37) 

Como consecuencia se obtiene la relacion de Josephson: 

a^2-vd. (1.38) 

Si bien la hipotesis de escala no ha sido demostrada con rigurosidad, esta 
surge de un niodo natural en el marco de las hipotesis de trabajo del grupo 
de renormalizacion. En este contexto vamos a describir la Ifnea argumental 
de Kadanoff que Ueva a la plausibilidad de la hipotesis de escala. 

Supongamos un Hamiltoniano tipo Ising como el descrito por la accion 



( 1.13 ), es decir un sistema de N espines en una red de dimension d. Con- 
sideremos la red dividida en celdas de lado aL, tenemos por tanto n = N/L'^ 
celdas cada una con L'' espines. Los argumentos de Kadanoff se aplican 
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solo a la situacion en que ^ 3> aL, es decir, en un entorno suficientemente 
pequeno del punto critico. En estas condiciones dentro de un misnio islote 
de espines correlacionados hay un gran niimero de celdas. 

Vamos a asociar a cada celda a (a = 1, 2, . . . , n) un espi'n CTq de acuerdo 
con la regla de la mayoria. Puesto que cada celda esta enteramente dentro 
de un islote de espines correlacionados, es razonable asumir que los cFa se 
comportaran como los espines individuales Ui en el sentido de que actuaran 
conio si tomasen los valores ±1. 

Podemos pensar que la accion escrita en terminos de estos momentos 
da , tenga la misma forma que la accion original, excepto que los parametros 
P y h seran en general distintos. Asociamos pues unos nuevos /3 y /i a la 
accion construida con las a^- Puesto que la temperatura cri'tica es una 
medida de la intensidad de la interaccion J, no haremos la discusion en 
funcion de /3 sino de la temperatura reducida e. 

Ahora bien, si la accion tiene la misma forma funcional, es razonable 
pensar que lo mismo ocurrira con los potenciales termodinamicos. Asi por 
ejemplo el potential de Gibbs de cada celda, G(e, h), sera la misma funcion 
de h y e que el potential de Gibbs correspondiente a cada site. Por ser este 
una magnitud extensiva se tiene: 

G{e,h) = L''G{e,h) . (1.39) 

La dependencia de h con h y e con e se puede pensar que sera en general 
una funcion de L: 

h = H{L)h , . 

l^T{L)e. ^ > 

En estas condiciones ya se puede argumentar el cumplimiento de la 
hipotesis de escala, sin embargo siguiendo con el argumento de Kadanoff, 
se asume que H(L) y T(L) se pueden expresar como una potencia de L: 
H{L) = y T{L) = , con x e y arbitrarios. Esto nos permite reescribir 



( 1.39 ) como: 

G(Lye,L''h) = L'^G{e,h) 

(1-41) 

=^G{Lyl'^e,L''/'^h) = LG{e,h) , 



que tiene la misma forma funcional que ( 1.34 ) excepto que L no es un 



parametro libre como lo es A en ( 1.34 ). En efecto L esta sujeto a la condition 
1 <C L^/a, sin embargo, puesto que ^ diverge en el punto de la transition 
sienipre podemos ponernos tan proximos al punto critico como queramos 
para lograr que L este en ese intervalo. Es decir que (|l.4l| ) significa que el 
potential de Gibbs es una funcion homogenea generalizada. 

Una Ifnea de argumentation analoga aplicada a la funcion de correlation 
a dos puntos, r(e,r), muestra que esta es tambien una funcion homogenea 
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generalizada. Teniendo en cuenta que la longitud de correlacon es funcion 
solo de la temperatura reducida, podemos escribir r(e,r) sustituycndo la 
dependencia en e por la dependencia en ^: r(^,r). Puesto que es una 
funcion homogenea tenemos: 

r(AC,Ar) =A"r(^,r) . (1.42) 

Tomando A = 1 /^: 

r(l,r/e) =r"r(C,r) . (1.43) 

Puesto que uno de los dos parametros esta fijo podemos escribir la 
funcion de correlacion como dependiente de un solo parametro, es decir 
r(l,r/^) = F{^,r). Es decir que podemos escribir: 

T{tr)=er{l,r/0 . (1.44) 

Esta igualdad significa que la correlacion entre dos espines situados a 
distancia r, depende de r solo a traves del cociente r/^. Dicho de otra 
forma, que solo existe una longitud caracterfstica en el problema y que esta 
es la longitud de correlacion. 

Sin embargo a la hora de extraer informacion mediante metodos numericos 
hay que tener en cuenta que las divergencias que marcan la existencia de 
una transicion de fase ocurren solo cuando el mimero de grados de libertad 
tiende a infinito. En un volumen finito, es decir dentro de los ordenadores, 
todas las cantidades se obtienen mediante sumas finitas siendo por tanto 
funciones analfticas en todo el espacio de acoplamientos. Como conse- 
cuencia aparecen dos fenomenos fundamentalmente: las divergencias son 
sustituidas por picos de altura finita, creciente con el tamafio de la red; la 
localizacion de estos picos se mueve en el espacio de acoplamientos conforme 
el tamafio de la red cambia. 

La teorfa de Finite Size Scaling (FSS) jlll esta basada en que estos 
efectos que aparecen al aumentar el tamafio de la red son precursores del 
comportamiento del sistema en el limite termodinamico, y que pueden ser 
explotados para extraer las propiedades de la transicion de fase. Considere- 
mos un retfculo de lado L. La teorfa de FSS descansa sobre la hipotesis de 
que en el entorno del punto crftico el comportamiento del sistema se puede 
describir enteramente en funcion de la variable reescalada: 

y^L/ar). (1.45) 

En una transicion de fase continua ^(T) crece conforme nos acercamos a 
Tc- Es razonable suponer que los efectos de tamafio finito se haran palpables 
cuando esta alcance el tamafio de la red L. Definimos la temperatura crftica 
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aparente T* (L) como aquella en la que £,{T* (L)) ^ L. En estas condiciones 
podemos estimar como depende de L la temperatura aparente puesto que: 



e(T;(i))(xi ^ |T;(i)-T,r 

=>\T*{L)~Tc\ - L-^/"" . (1.46) 

En general si Ql{T) es un observable fi'sico en el sistema finito, cuyo 
comportamiento en el limite termodinamico es: 

Qoo{T) ^ Coo^'P e^O. (1.47) 

En un volumen finito el comportamiento sera del tipo: 

Ql{T) ^ CLe-P , (1.48) 

cone= \T~T*{L)\ - L^i/'". 

Si introducimos ahora la hipotesis de FSS se tiene: 

QL{T)^L^f{y) . (1.49) 



Para determinar uj basta imponer que ( 1.49| ) reproduzca (1.47) en el 
limite de L — !■ cx): 

f{y) - Cooe-P i -> oo , (1.50) 

para ello u) = p/v. 

Por otra parte puesto que en el sistema finito no hay transicion de fase: 

/(2/)->/o e^O, (1.51) 

es decir f[y) debe aproximar una constante finita en el punto cn'tico. 

Para resumir, el resultado significativo de este analisis, desde el punto 
de vista de extraer informacion numerica, es que la forma en que varfan 
las magnitudes termodinamicas con el tamafio de la red esta determinada 
por los exponentes cn'ticos del sistema en el limite termodinamico: dado un 
observable Q{e) que en L = oo diverge en el punto de la transicion como 
(3(e) ^ e-^P , en L finito su comportamiento sera de la forma: 

Ql{T:{L))^ f^LP''^ . (1.52) 

1.3.2 Transiciones de fase de primer orden. 

Como ya se ha sefialado anteriormente, segiin el criterio basico de clasifi- 
cacion termodinamica las transiciones de fase de primer orden son aquellas 
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en las que la primera derivada de la energia libre presenta una discon- 
tinuidad. 

Cuando un sistema realiza una transicion de fase de primer orden de 
una fase de alta temperatura a otra de baja temperatura, este desprende 
una cantidad no nula de calor (el calor latente C;) a medida que se enfria 
en un intervalo infinitesimal de temperaturas alrededor de la temperatura 
de la transicion. El valor de C/ se relaciona con el salto en la entropia que 
se produce en las transiciones de primer orden de la siguiente forma: 

n^-{dF/dT) =^ Ci=T,An. (1.53) 

Se tiene mucho menor conocimiento teorico sobre el comportamiento 
de las transiciones de fase de primer orden que sobre los fenomcnos aso- 
ciados a las transiciones continuas. Esto es debido a que al no haber una 
longitud de correlacion divergente, no podemos restringir el estudio a los 
fenomenos asociados a largas longitudes de onda, por lo tanto el concepto 
de Universalidad no es aplicable. Las singularidades que aparecen en las 
transiciones de primer orden se deben exclusivamente a la coexistencia de 
fases, no hay region crftica, ni exponentes cri'ticos. No podemos usar por 
tanto la teorfa de FSS para estudiar el comportamiento de estas transiciones 
en un volumen finito. 

Debido a la existencia de calor latente, el calor especi'fico es una funcion 
5{e). Analogamente a lo que ocurre con las transiciones continuas, en una 
red finita de lado L, las singularidades de las funciones en el punto de la 
transicion, en este caso las (5's, seran suavizadas y sustituidas por picos 
de altura finita. De la misma forma esperamos un desplazamiento de la 
temperatura a la cual aparecen estos picos con el tamaiio de la red, T*[L). 

Para discutir cuantitativamente los efectos de tamafio finito asocia- 
dos a una transicion de primer orden vamos a utilizar la aproximacion 
fenomenologica introducida por Binder [ p^ . Se asumira que el tamaiio de 
la red es mayor que la longitud de correlacion de cada fase, de tal manera 
que los observables toman los valores asintoticos de la transicion. 

En una red de lado L y a, una cierta temperatura T, la distribucion de 
probabilidad de la energia, Pl{E), es gausiana: 

(1.54) 

donde Eq es la energia del sistema en el limite de volumen infinito, y 
es caracterfstica de la temperatura T. La anchura de la distribucion es 
proporcional al calor especi'fico en este Ifmite, denotado por C. 

En el punto de la transicion de fase la caracterfstica fundamental de las 
transiciones de primer orden es la coexistencia entre fases. Supongamos 
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que ocurre una transicion enT = Tc de un estado desordenado a un estado 
ordenado con degeneracion q. Denotemos por (i?-) la energia interna 
de la fase desordenada (ordenada) aT = Tc, de manera que el calor latente 

es E, - E^. 



La gencralizacion de ( 1.54 ) al caso en que hay coexistencia de fases es 
considerar que la distribucion de probabilidad a Tc es una doble gaussiana. 
En el intervalo AT — T — Tc las gaussianas estaran centradas respectiva- 
mente en E+ + C+AT y E^ + C-AT: 



P (p^ / {E-E+-C+AT)^L\ 
Pl{E) oc ^==cxp{ 1 (1.55) 



, (E ^ E^ -C-ATfL'^ ^ 
exp{-^ } , (1-56) 



siendo C+ y C_ el calor especffico de las fases desordenada y ordenada 
respectivamente: 

lim a(T) = C_ (1.57) 
lim C„(T) = C+ , (1.58) 

Puesto que se esta considerando el comportamiento del sistema cerca 
de Tc se puede considerar que C+ y C_ son constantes a lo largo de esta 
region. Los dos picos estan pesados de acuerdo con la diferencia de energia 
libre entre las dos fases AF{T) = F+ ^ F^, de manera que lejos de Tc 
solo uno de los picos sobrevive, mientras que en el punto de la transicion 
AF(Tc) = y ambos terminos contribuyen con el mismo peso: 

(1.59) 
(1.60) 




donde 



(1.61) 



2T 



Se puede aproximar x expandiendo en serie F^ y F^ entorno a Tc hasta 
orden AT (orden mayor significarfa considerar la variacion de C+ y C_). 
Teniendo en cuenta que F± = U± — TS± y que dF± = —S±dT se obtiene: 

AT ^ -(^^-^-)^^ (1.62) 
Tc ^ ' 
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El valor medio de la energia vendra dado por el primer momento de la 



distribucion de probabilidad (1.56) 



(£;) = "+^+ + "-^-+Ar°+^+ + "-^' , (1.63) 

y por lo tanto el calor espccifico: 

Cv[T,L) = = — — (1.64) 



dT a+ + a 

a+a^L'^ [{E+ - E^) + {C+ - C - ) AT] 
^ (a+ + a-)2 

cuyo maximo ocurre para: 

T*{L)-Tc ^ T.logg 1 

con mia altura de: 



2 



(1.65) 



(1.66) 



a(i)-- = ''^^^^L^ + • (1-67) 

Cuando se esta interesado en simulaciones de Monte Carlo hay errores 
asociados a la simulacion que pueden hacer mas aceptable una descripcion 
de los efectos de tamano finito en funcion de exponentes efectivos: 

T*{L)-Tc{oo) cx L-^ (1.68) 
C;°^'^(L) (X i"" , (1.69) 

dondc hemos introducido los indices A y am para estudiar esta depen- 
dencia con L. Es decir, el resultado que se encuentra es que los exponentes 
A y am son igual a la dimension del espacio, d, asintoticamente, es decir 
para L suficientemente grande (L ^ ^). 

Sin embargo la definicion de estos exponentes es puramente formal, y no 
tienen nada que ver con las definiciones que hemos visto En una transicion 
de primer orden L aparece solo porque el volumen es L'^ en d dimensiones: 
los maximos de las susceptibilidades y del calor especffico crecen como L'^ 
y la funcion 5 se obtiene en el h'mite de volumen infinito porque la anchura 
de estas funciones decrece como L~'^. 

Desde el punto de vista del estudio numerico, en las transiciones de 
primer orden, en general, cualquier observable termodinamico muestra en 
la evolucion de MC saltos entre las dos fases {Bip-flops) generando una 
distribucion de doble pico. 
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Las transiciones fuertes (pequena ^) son faciles de detectar. Las discon- 
tinuidades se observan ya a nivel de los ciclos de histeresis que muestran ra- 
mas metaestables. Las discontinuidades en la magnitudes termodinamicas 
son relativamente faciles de detectar, y tanto el calor especifico como las 
susceptibilidades en transiciones magneticas divergen como para redes 
de tamaiio no muy grande. 

La situacion numcrica cs menos clara en las Uamadas transiciones do 
primer orden debiles. La longitud de correlacion es muy grande, y por tanto 
los efectos transitorios se prolongan hasta redes de tamaiio respetable para 
nuestro actuales ordenadores. En estos casos uno dcbe simular redes de 
L creciente y tratar de ver si el comportamiento de primer orden anteri- 
ormente descrito se alcanza en el Hmite de volumenes muy grandes. En 
particular se pueden usar las tecnicas desarroUadas de FSS para estudiar 
los exponentes crfticos efectivos y ver como evolucionan con L las magni- 
tudes termodinamicas. Dentro de este contexto se dice que en la region 
asintotica las transiciones de primer orden tienen asociado un exponente 
critico "z/"= 1/do que "a/z^" = d, pero entendiendo bien que no son expo- 
nentes crfticos definibles en sentido extricto. 
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Chapter 2 

El Modelo 0(4) 
Anti-Ferromagnetico 
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2.1 Introduccion 



En cl Lagrangiano del Modelo Estandar el campo de Higgs esta descrito por 
un campo escalar de cuatro componentes, el cual, aparte de interaccionar 
con los campos gauge y con los fermiones, interacciona consigo mismo a 
traves de un autoacoplo cuartico. Si consideramos el Ifmite en el que los 
acoplos gauge y los de Yukawa tienden a cero, los campos gauge y los 
fermionicos se desacoplan. La accion correspondiente a este li'mite es el 
modelo A*^: 

5eo„t = [ I 5](a^*B(r))^ + im|*B(r)2 + f ^bW^ . (2.1) 

Una forma conveniente de escribir la accion cuando estamos interesados 
en simulaciones de MC o expansiones de alta T es: 

^lat-^(-«^*r*r+/i +*r+A[*?-l]') , (2.2) 
r ^l 

donde n recibe el nombre de hopping parameter. 

La correspondencia entre los parametros en la red y en el contfnuo es la 
siguiente: 

*B(r)a = V2K*r ; 

TO|a2 = (1-2A)/k-8; (2.3) 
gB = 6A/k2 . 

El diagrama de fases del modelo ( |2.2[ ) presenta una linea de transiciones 
de segundo orden Kc(A) separando la fase simetrica de la fase en que la 
simetrfa 0(4) esta rota a 0(3)(8)0(3). El li'mite contfnuo hay que tomarlo 
aproximandose a esta li'nea crftica. 

En A = se tiene el punto fijo Gaussiano infrarojo puesto que la teorfa 
en este limite se reduce al campo escalar libre. El Hmite contfnuo en el 
entorno de este punto fijo es trivial en el sentido de que la constante de 
acoplamiento renormalizada se hace cero. 

Para poder definir un Ifmite contfnuo no trivial (gR ^ 0) es necesario 
encontrar un punto fijo ultravioleta. 

En dimension d > A esta probado rigurosamente ||l|, D que A^'' no tiene 
ningiin punto fijo ultravioleta, la teorfa es por tanto trivial para todo valor 
de A. 

En dimension d = 4 las evidencias analfticas ^, §| y numericas (|, |6[ 
acumuladas soportan la conjetura de Wilson apuntando a la trivialidad 
de las teorfas A**. Sin embargo una prueba rigurosa como la obtenida en 
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c? > 4 no ha sido obtenida, a pesar de los serios intentos Uevados a cabo 
(ver por ejemplo Q). 

Los teoremas existentes cubren solo ciertos casos especiales, que por ar- 
gumentos de plausibilidad, se supone que las QFT deben cumplir En 
particular se asume generalmente que la autointeraccion del campo escalar 
se puede tratar perturbativamente. Sin embargo las teorfas con acoplos 
negativos generan fuertes oscilaciones en los campos, que podrfan dar lugar 
a que se generen no perturbativamente nuevos puntos fijos, donde even- 
tualmente se podrfan definir Ifmites contmuos no triviales [^] . Desde este 
punto de vista el estudio de teorfas Anti-ferromagneticas (AF) es intere- 
sante puesto que nada impide formular una una QFT, "honesta" en el sen- 
tido de que cumple los axiomas de Osterwalder-Schrader, usando modelos 
antiferromagneticos (AF) (ll]. 

El antiferromagnetismo ha sido considerado en una amplia variedad de 
modelos con el objetivo de encontrar propiedades no presentes en teorfas 
puramente ferromagneticos. En el contexto de la superconductividad de 
alta temperatura el AF parece jugar un papel esencial |p^ . En dimension 
d = 4, interacciones competitivas como posibles agentes de nuevas clases 
de universalidad han sido investigadas para estudiar el punto multicrftico 
de modelos tipo Yukawa. Este punto multicrftico es el punto de encuentro 
de cuatro fases distintas (FM, AF, Ferrimagnetica y Paramagnetica (PM)). 
La cuestion de si es posible o no definir un Ifmitc contfnuo no trivial en este 
punto permanece como un problema abierto ||l3|, |l^ . 

Asf pues el estudio de modelos AF se presenta como una posibilidad 
interesante de estudiar en detalle. Estos modelos suelen presentar diagra- 
mas de fases muy ricos, y presumiblemente nuevas clases de universalidad 
podrfan ser encontradas JT?! . La existencia de acoplamientos de signo o- 
puesto influencia el vacfo de la teorfa, en concreto en el estado fundamental 
se pueden observar fenomenos de frustracion (la energfa no puede ser mi- 
nimizada simultaneamente para todos los acoplamientos) o de desorden (la 
entropfa del vacfo es distinta de cero). 

En el trabajo que se presenta a continuacion se considera la teorfa 
en el If mite en que A — > cxd, lo cual equivale a fijar el modulo del campo 
de Higgs. Para introducir Antiferromagnetismo se ha afiadido un acoplo 
negativo a segundos vecinos. El objetivo es la biisqueda y caracterizacion 
de puntos de transicion de fase de segundo orden en el diagrama de fases. 
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Abstract 



We study the phase diagram of the four dimensional 0(4) model with 
first (Pi ) and second ((32 ) neighbor couplings, specially in the (32 < region, 
where we find a line of transitions which is compatible with second order. 
We also compute the critical exponents on this line at the point /?i = (F4 
lattice) by Finite Size Scaling techniques up to a lattice size L — 24, being 
these exponents different from the Mean Field ones. 

2.2 Description of the model 

Our starting point is the non-linear a model, with action: 

Sa ^ -(3Y,^r^r+il ■ (2.4) 

Where $ is a 4-coniponent vector with fixed modulus 4>r • 4>r = 1- 
The naive way to introduce AF in the non-linear a model is to consider 
a negative coupling. In this case the state with minimal energy for large /3 is 
a staggered vacuum. On a hypercubic lattice, if we denote the coordinates 
of site r as {rx,ry,rz,rt), making the transformation 

*r ^ (-!)'■-+'■«+'■-+'■* *r , (2.5) 

the system with negative /3 is mapped onto the positive P one, both regions 
being exactly equivalent. 

Therefore to consider true AF we must take into account either different 
geometries or more couplings, in order to break the symmetry under the 
transformation ( |2.5| ). In four dimensions the simplest option is to add more 
couplings, we have chosen to add a coupling between points at a distance 
of V2 lattice units. 

Following this we will consider a system of spins {$r} taking values in 
the hyper-sphere S'^ C and placed in the nodes of a cubic lattice. The 
interaction is defined by the action 

5 = ^ *r*r+A - /32 ^ *r*r+A+P , (2-6) 

The transformation ( ^.5[ ) maps the semi-plane (3i > onto the (3i < 0, 
and therefore only the region with (3i > will be considered. On the line 
(3i = the system decouples in two F4 independent sublattices. 

When (32 = the model is known to present a continuous transition 
between a disordered phase, where 0(4) symmetry is exact, to an ordered 
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phase where the 0(4) symmetry is spontaneously broken to 0(3). This 
transition is second order, being the critical exponents those of MFT: a = 0, 
V = 0.5, /3 = 0.5, 77 = and 7 = 1 up to logarithmic corrections. The critical 
coupling for this case can be studied analytically by an expansion in powers 
of the coordination number (5 = 2d), being ji'^ = 0.6055 + 0{q~'^'^) |T^ . 

From a Mean Field analysis, we observe that for P2 > the behavior 
of the system will not change qualitatively from the (^2=0 case but with 
higher coordination number. In fact, taking into account that the energy 
(for non-frustrated systems) is approximately proportional to the coordi- 
nation number, there will be a transition phase line whose approximate 
equation is 

/3J + Q(3-2 = , (2.7) 

where Q is the quotient between the number of second and first neighbors, 
2d{d — 1) and 2d, respectively. This line can be thought as a prolongation 
of the critical point at /32 = so the transitions on this line are expected 
to be second order with MFT exponents. This is also the behavior of the 
two couplings Ising model in this region 

When (32 < 0, the presence of two couplings with opposite sign makes 
frustration to appear, and very different vacua are possible. 



2.3 Observables and order parameters 

We define the energy associated to each coupling: 

^l-^-E*-*r+,, (2.8) 



- ^ - E . (2.9) 



In terms of these energies, the action reads 

S = -p,Ei - P2E2 . (2.10) 
It is useful to define the energies per bound as 

ei = — El, 62 = E2 , (2.11 

AV 12V 

where V = L'^ is the lattice volume. With this normalization ei , 62 belong 
to the interval [—1,1]. 
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We have computed the configurations which minimize the energy for 
several asymptotic values of the parameters. We have only considered con- 
figurations with periodicity two. More complex structures have not been 
observed in our simulations. 

Considering only the /3i > case, we have found the following regions: 

1. Paramagnetic (PM) phase or disordered phase, for small absolute 
values of /?2- 

2. Ferromagnetic (FM) phase. It appears when (3i + 6/32 is large and 
positive. 

When the fluctuations go to zero, the vacuum takes the form = v, 
where v is an arbitrary element of the hypcr-sphere. 

Concerning the definition of the order parameter let us remark that 
because of tunneling phenomena in finite lattice we are forced to use 
pseudo-order parameters for practical purposes. Such quantities be- 
have as true order parameters only in the thermodynamical limit. In 
the FM phase, we define the standard (normalized) magnetization as 

r 

and we use as pseudo-order parameter the square root of the norm of 
the magnetization vector 

Mf = (^Ml) . (2.13) 

This quantity has the drawback of being non-zero in the symmetric 
phase but it presents corrections to the bulk behavior order 1 / W. 

3. Hyper-Plane Antiferromagnetic phase (HPAF). It corresponds to large 

/?!, with (32 in a narrow interval ([— /3i/2, — /3i/6] in the Mean Field ap- 
proximation). In this region the vacuum correspond to spins aligned 
in three directions but anti-aligned in the fourth (/i). 

In absence of fluctuations the associated vacuum would be $r = 

{—lyt^v, where ^ can be any direction, and v any vector on S'*. 

We define an ad hoc order parameter for this phase as 

Mhpaf,^ = ^ Y.^-lY'-^r . (2.14) 

r 

Mhpaf,/* will be different from zero only in the HPAF phase, where 
the system becomes antiferromagnetic on the /x direction. Prom the 
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four order parameters (one for every possible value of /i) only one of 
them will be different from zero in the HPAF phase. So, we define as 
the pseudo order parameter: 



AWf-^^MIp^p^^ . (2.15) 

4. Plane Anti-Ferromagnetic (PAF) phase for /32 large and negative. In 
this region the ground state is a configuration with spins aligned in 
two directions and anti-aligned in the remaining two. It is character- 
ized with by one of the six combinations of two different directions 
(/i,i^), and an arbitrary spin v: $r = (— l)'"'' ^'^'^ v. For the PAF 
region we first define 

MpAF,p.. = ^ ^(-l)'--+'^"*r , (2.16) 



and the quantity we measure is 



MpAF- /Em|af.(,,.) (2-17) 

In order to avoid undesirable (frustrating) boundary effects for ordered 
phases, we work with even lattice side L as periodic boundary conditions 
are imposed. 

From this data we can compute the derivatives of any observable with 
respect to the couplings as the connected correlation function with the 
energies 

do 



{OE,)-{0){E,) (2.18) 



An efficient method to determine f3c for a second order transition is to 
measure the Binder cumulant pO| for various lattice size and to locate the 
cross point in the space of /3. 

For O(A^) models Ul{P) takes the form Q: 



C/i(/3) = l + 2/iV-A^-f/ (2.19) 



where m is an order parameter for the transition. 

It can be shown [|o[ |2l) that Ul{0) 0{1/V) and Ul{oo) 2/N. 
The slope of Ul{(3) at /3c increases with L. 
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The value of the Binder cumulant is closely related with the triviality 
of the theory since the renormalized coupling (in the massless thermody- 
namical limit) at zero momentum can be written as: 

gn = lim gR{L) = lim {L/^LfUL{f3c) (2.20) 

L — >oo L — >oc 

where is the correlation length in the size L lattice. 

From this point of view triviality is equivalent to have a vanishing 
in the thermodynamical limit. In this context it is clear that we can use 
the value of g^ to classify the universality class. Out of the upper critical 
dimension, i/^L is a constant at /3c since ^ ~ i, and we could use the 
Binder cumulant for the same purpose p2| . At the upper critical dimension, 
£,L presents logarithmic corrections and L/^l is no longer a constant at 
Pc- For the FM 0(4) model in d = 4 (upper critical dimension) we have 
perturbatively L/S,l ~ (Ini)"^/"* |Q. In order to have a non trivial theory, 
the Binder cumulant should behave as a positive power of In L, but from its 
definition we see that Ul{P) < 1- This is just another way of stating 
the perturbative triviality of the FM 0(4) model. 



2.3.1 Symmetries on the F4 lattice 

In the f3i — case the system decouples in two independent lattices, each 
one constituted by the first neighbors of the other. So we consider two 
lattices with F4 geometry. There are several reasons to choose the point 
Pi = for a careful study of the PM-PAF transition. The region with 
Pi > 1.5 evolve painfully with our local algorithms; For small Pi we expect 
very large correlation in MC time because the interaction between both 
sublattices is very small, and the response of one lattice to changes in 
the other is very slow. We also remark that the presence of two almost 
decoupled lattices is rather unphysical. 

We also have the experience from a previous work for the Ising model 
ijl^ that the correlation length at its first order transition is smaller in the 
F4 lattice, that means, we can find asymptotic critical behavior in smaller 
lattices. 

However we should point out that the results in the F4 lattice cannot 
be easily extrapolated to a neighborhood of the Pi axis. Certainly, the 
geometry of the model is very modified when pi 0, and perhaps continuity 
arguments present problems. Nevertheless, we have run also the case Pi ^ 
0, and as occurs in the Ising model we have not found qualitative differences. 

In the following when we refer to the size of the lattice L on the F4 
lattice we mean a lattice with sites. 

We have to find the configurations that maximize E2 in order to define 
appropriate order parameters for the phase transition. 
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The system has a very complex structure. As starting point we have 
studied numerically the vacuum with /32 -C 0. For this values we have found 
in the simulation: 

1. The vacuum has periodicity two. To check this, we have defined: 

Vi = Z^E*^- (2-21) 

where i = 0, . . . , 7 stands for the i*'' vertex of each 2'* hypercube 
belonging to the F4 lattice, and with I we denote the 2* hypercubes 
themselves. 

From these vectors we can define the 8 magnetizations associated to 
the elementary cell, 

Vi = (/V^) , (2.22) 

We have checked that all Vt tends to 1 for the ordered phase in the 
thermodynamical limit, so we conclude that the ordered vacua have 
periodicity two. 

Let us remark for the sake of completeness that all order parameters 
we have defined can be written as an appropriate linear combination 
of the Vi. 

2. In the elementary cell, ^r+p+P = V/i, ly with ji < v. So, in this 
section we will restrict the study of the vacuum structure to the four 
sites {i = 0, 1, 2, 3) belonging to the cube in the hyper-plane rt = 0. 

3. We have measured the energy per bound associated to the second 
neighbors coupling. We check that in the thermodynamical limit 

62 = -1/3. 

4. If we choose the symmetry breaking direction by keeping fix one vec- 
tor, (eg. $0) we find: 

3 

5^ ((*o • *i)*o - *i) = , (2.23) 

i=l 

The vacuum structure is not completely fixed by these three conditions 
since different symmetry breaking patterns are possible. For instance, a 
configuration *o = (1,0,0,0), *i = (-1/3, ^Vi), *o = (-1/3, ^Va), 
$0 = (^1/3, ^^va), with Vi a 3-component unitary vector with the con- 
straint X^j^j ViVj = 0, breaks 0(4), but an 0(2) symmetry remains (for 
the different Vi). 
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To determine which is the vacuum in presence of fluctuations, we con- 
sider four independent fields in a 2^ cell with periodic boundary conditions. 
Let us first consider an 0(2) group. We can study the four vectors as a me- 
chanical system of mass- less links of length unity, rotating in a plane around 
the same point, whose extremes are attached with a spring of natural length 
zero. The energy for the system is: 



E = - ^ cos{9i - Oj) . 

ij=0,i>j 



(2.24) 



We consider the fluctuation matrix, H = dE"^ /dOidOj in order to find 
the normal modes. The matrix elements of H take the form: 

= Sij J2 cos(^i - Ok) - cosiOi - ej){l - Sij) , (2.25) 

In the FM case the minimum correspond to 6i = <j), for all i. The 
fluctuation matrix is 



FM 



/ -3 1 1 1 \ 

I- 311 

II- 31 
VI 1 1-3/ 

which has a single zero mode, and a three times degenerated non-zero mode 
with eigenvalue A = — 4. 

For the AF (maximum energy) case, the maximum energy is found, up 
to permutations, dX 9q = (j), 9i = (j) + ir,62 = (f' + a and 6*3 — (fi + n + a, Va. 
In addition to the (f> freedom that corresponds to the global 0(2) symmetry, 
there is a degeneration of the vacuum in the a angle and this zero mode is 
double Va. 

The fluctuation matrix in the AF case is 



Haf = 



( 1 -1 

-1 1 

cos a — cos a 

\ — cos a cos a 



cos a — cos a \ 

— cos a cos a 

1 -1 

-1 1 



The other two eigenvalues are Ai,2 = 2(1 ± cos a), so, an additional zero 
mode appears when a = 0, obtaining in this case a three fold degenerated 
zero mode corresponding to: 60 = Oi = 62 + tt = O3 + w. 

The 0(4) case is qualitatively similar. Wc have 12 degrees of freedom. 
Of all configurations that minimize the energy, that with a largest degener- 
ation (9-times) consist of 2 spins aligned and 2 anti-aligned that correspond 
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to a PAF vacuum. We consider this degeneration as the main difference 
with the FM sector, and could be relevant to obtain different critical expo- 
nents. 

In presence of fluctuations the configurations with largest degeneration 
are favored by phase space considerations, so we expect that the real vac- 
uum is a PAF one. This statement will be checked below with Monte Carlo 
data in the critical region. 



2.4 Finite Size Scaling analysis 

Our measures of critical exponents are based on the FSS ansatz ^ . 
Let be {0{L,P)) the mean value of an observable measured on a size L 
lattice at a coupHng /3. If O(oo, (3) ~ |/3 — (3c\^° , from the FSS ansatz one 
readily obtains 

(0(L, /?)) = L-^'/'-FoiL/Ci^, /?)) + ..., (2.26) 

where Fq is a smooth function and the dots stand for corrections to scaling 
terms. 



To obtain u we apply equation (2.26) to the operator dlogMpAF/d/3 



whose related x exponent is 1. As this operator is almost constant in the 
critical region, we just measure at the extrapolated critical point or any 
definition of the apparent critical point in a finite lattice, the difference 
being small corrections-to-scaling terms. 

For the magnetic critical exponents the situation is more involved as the 
slope of the magnetization or the unconnected susceptibility is very large 
at the critical point. 

We proceed as follows (see refs. |2^ for other applications of this 
method). Let be any operator with scaling law xq ~ 1 (for instance 
the Binder parameter or a correlation length defined in a finite lattice di- 
vided by L). Applying eq. ( ^.26 ) to an arbitrary operator, O, and to Q we 



can write 

(0(i, /3)) = L-«/Vo,0((e(L, /?))) + .... (2.27) 

Measuring the operator O in a pair of lattices of sizes L and sL at a coupling 
where the mean value of Q is the same, one readily obtains 



(0(sL,/3)) 



(0(L,/3)) 



= s^o/^-f... . (2.28) 



The use of the spectral density method (SDM) ||2^] avoids an exact a priori 
knowledge of the coupling where the mean values of 8 cross. We remark 
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that usually the main source of statistical error in the measures of magnetic 
exponents is the error in the determination of the coupling where to mea- 



sure. However, using eq. (2.28) we can take into account the correlation 
between the measures of the observable and the measure of the coupling 
where the cross occurs. This allows to reduce the statistical error in an 
order of magnitude. 

2.4.1 FSS at the upper critical dimension: logarithmic 
corrections 

Being d = 4 the upper critical dimension of the FM 0(4) model, logarith- 
mic corrections to the Mean Field predictions are expected to set in. In 
particular, FSS in its standard formulation breaks at c? = 4 because the 
essential assumption, namely S,l{Pc) i is no longer true. In fact, in four 
dimensions 

^{oo,t) ^ \t\-^/''\ln\t\\i . (2.29) 

The FSS formula for the correlation length was calculated by Brezin 
|^3|] . At the critical point one gets: 

C(L,/?e)^i(lni)i/^ (2.30) 

It has been suggested j28| that the usual FSS statement should be re- 
placed by the more general one: 

When applying the quotient method described above to systems in four 
dimensions one has to take into account the logarithmic corrections, so that 
the modified formula reads: 



{0{sL,l3)) 



{0{L,P)) 



e(L,f3)=eisL,l3) 



1 \ 1/4 

s^^/'^fl + j^j . (2.32) 



This point is particularly important when measuring the magnetic crit- 
ical exponents because as we have already mentioned, the slope of the 
magnetization and susceptibility are very large, and special care has to be 
taken when locating the coupling where to measure. 

2.5 Numerical Method 

We have simulated the model in a L'* lattice with periodic boundary con- 
ditions. The biggest lattice size has been L = 24. For the update we have 
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employed a combination of Heat-Bath and Over-relaxation algorithms (10 
Overrelax sweeps followed by a Heat-Bath sweep) . 

The dynamic exponent z we obtain is near 1. Cluster type algorithms 
are not expected to improve this z value. In systems with competing in- 
teractions the cluster size average is a great fraction of the whole system, 
loosing the efficacy they show for ferromagnetic spin systems. 

We have used for the simulations ALPHA processor based machines. 
The total computer time employed has been the equivalent of two years of 
ALPHA AXP3000. We measure every 10 sweeps and store the individual 
measures to extrapolate in a neighborhood of the simulation coupling by 
using the SDM. 

In the F4 case, we have run about 2 x lO^r for each lattice size, being 
T the largest integrated autocorrelation time measured, that corresponds 
to MpAF, and ranges from 2.3 measures for L = 6 to 8.9 for L = 24. We 
have discarded more than lO^r for thermalization. The errors have been 
estimated with the jack-knife method. 

2.6 Results and Measures 
2.6.1 Phase Diagram 

We have studied the phase diagram of the model using a L = 8 lattice. 
We have done a sweep along the parameter space of several thousands 
of iterations, finding the transition lines shown in Figure 1. The symbols 
represent the coupling values where a peak in the order parameter derivative 
appears. 

The line FM-PM has a clear second order behavior. It contains the 
critical point for the 0(4) model with first neighbor couplings (/3i « 0.6, 
/32 = 0) with classical exponents {y = 0.5, r] = 0). In the /?i = axis, we 
have computed the critical coupling (/J^ ~ 0.18) and the critical exponents 
as a test for the method in the F4 lattice. We have also considered the 
influence of the logarithmic corrections when computing the exponents. 

The lines FM-HPAF, HPAF-PAF and PM-HPAF show clear metasta- 
bility, indicating a first order transition. 

The regions between the lower dotted line and the PAF transition line, 
and between the upper dotted line and the FM transition line, are disor- 
dered up to our numerical precision. We could expect always a PM region 
separating the different ordered phases, however, from a MC simulation it is 
not possible to give a conclusive answer since the width of the hypothetical 
PM region decreases when increasing /3i, and for a fixed lattice size there 
is a practical limit in the precision of the measures of critical values. 
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Figure 2.1: Phase diagram obtained from tlie MC simulation on a L = 8 lattice 



The line PM-PAF has a very interesting behavior. To get ride of the 
influence of the HPAF region we have started the study at /3i = 1.5. The 
energy distributions encountered at /3i = 1.5 and Pi = 1.0 are displayed in 



figure 2.2. We do not obeserve latent heat anymore in lattices up to L=16 



in /3i =0. Two scenarios are suggested by this fact: 

1. There exist a value of f3i in which the order of the phase transition 
changes from being first order, to be continuous. 

2. The phase transition line is first order everywhere, though increasily 
weak as the limit /3i = is approached. 

The possibility of a second order behavior of the PM-PAF transition line 
contrast with the first order one found in the Ising model with two couplings 



in the analogous region |19|. This would not be that surprising because 
we are dealing now with a global continuous symmetry. The spontaneous 
symmetry breaking of such symmetries manifest in the appearance of soft 
modes or low energy excitations (long wavelength), the Goldstone bosons 
in QFT terminology |p9|| . In general, these low energy modes will perturb 
the mechanism of long distance ordering, softening in this way the phase 
transitions. 
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Figure 2.2: Energy distribution for L=8,12,16 measured at the peak of the 
specific heat at f3i — 1.5 (upper window) and at /3i — 1.0 (lower window). 



Regarding the differences with the FM case, the most remarkable fea- 
ture is the different vacuum structures appearing, specially the very large 
degeneration in the FAF transition, in contrast with the single degeneration 
of the FM 0(4) mode. 

For the reasons mentioned in section 1.3.1, the simpler point for studying 
the properties of the transition, namely the critical exponents is the F4 
limit. Most of the MC work has been done for this case. 



2.6.2 Results on the F4 lattice 
Results on the FM region 

Firstly, we have checked our method on the FM region of the F4 lattice. 
In Figure 3 the crossing points of the Binder cumulant for various lattice 
sizes are displayed. The prediction for the critical coupling /3c ~ 0.1831(1) 
agrees with an earlier study by Bhanot |3^ . 

Concerning the measures of critical exponents, we have applied the quo- 
tient method, described in section 1.4 . In table |2.l| we q uote the results 
when logarithmic corrections are included (formula ( 2.32| )), and also for 
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sake of comparison, when they are neglected (formula (2.28)). We see how 
in fact the agreement of the critical exponents with the MFT predictions 
is better when the logarithmic corrections are taken into account. 



L values 


8/16 


12/16 


10/12 


(without logarithmic corrections) 


ajv 


0.08(5) 


0.02(2) 


0.13(12) 


[3/y 


0.92(3) 


0.94(3) 


0.87(4) 




2.16(2) 


2.12(2) 


2.24(4) 


(with logarithmic corrections) 


ajv 


0.0 


0.0 


0.03(8) 




1.04(3) 


1.06(3) 


1.04(2) 




1.94(3) 


1.90(4) 


1.93(3) 



Table 2.1: Critical exponents for the FM-PM phase transition in the F4 
lattice. 




-0.1125 -0.1100 -0.1075 -0.1050 -0.1025 -0.1000 

Figure 2.3: Energy distribution for L=16,20 and 24 at the peak of the specific 
heat on the F4 lattice. 
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Figure 2.4: Crossing points of the Binder Cumulant for various lattice sizes on 
the FM-PM phase transition 



From now on we will focus on the transition between the PM phase and 
the PAF phase on the F4 lattice. 

Vacuum symmetries on the PAF region 

We will check using MC data that the ordered vacuum in the critical region 
is of type PAF. 
Let us define 

Ay = V, • V, . (2.33) 

The leading ordering corresponds to the eigenvector associated to the max- 
imum eigenvalue of the matrix that should scale as L~'^^/" at the critical 
point. The scalin g la w of the biggest eigenvalue agrees with the (3/v value 



reported in Table 2.3, and the associated eigenvector is, within errors, (1,1, 
1,-1). 

We also have found that the other eigenvalues scale as L^"^. This is the 
expected behavior if just the 0(4) symmetry is broken, and it remains an 
0(3) symmetry in the subspace orthogonal to the 0(4) breaking direction. 
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Figure 2.5: Crossing points of tlie Binder cumulant for various lattice sizes. 



Critical Coupling 

To obtain a precise determination of th e critical point, f3c, we have used 
the data for the Binder parameter ( ^.19| ). 

In Figure 4 we plot the crossing points of the Binder cumulants for the 
simulated lattices sizes. Extrapolations have been done using SDM from 
simulations at {32 = -0.7090 for L = 6,8, 10, 12 and 16; P2 = -0.7078 for 
L = 20, and 02 = -0.7070 for L = 24. 

The shift of the crossing point of the curves can be explained through 
the finite-size confluent corrections. The dependence in the deviation of the 



crossing point for L and sL size lattices was estimated by Binder 1 20 



/3c(L,sL)-/3c 



1 



^ — CiJ— 1 / u 



where to is the universal exponent for the corrections-to-scaling. 
The infinite volume critical point the value 

/3c = -0.7065(5) [-f 2] [-2] , 



(2.34) 



(2.35) 



where the errors in brackets correspond to the variations in the extrapola- 
tion when we use the values uj — 0.5 and w = 2 respectively. In Figure 5 
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Figure 2.6: Extrapolation to /3c(oo) for Li — 6, 8, 10, 12 (circles, crosses, triangles 
and squares symbol respectively) 



we plot eq. ( p. 34 ) for u = 1 



Using the previous value of (3c we can compute the Binder cumulant at 



this point. In table 2.2 we quote the obtained values. The result points to 
that the Binder cumulant stays constant in the critical region. This result 
would be compatible with a non zero value of the renormalized coupling 
when L increases. 

Concerning the possibility of having logarithmic corrections in the de- 
termination of the critical coupling, from the numerical point of view, it is 
not possible to discern between the to effect, and a logarithmic correction. 

Thermal Critical exponents: a, v 

The critical exponent associated to correlation length can be obtained from 
the scaling of: 

where M is an order parameter for the transition, AfpAF for our purposes. 
In the critical region n ~ L^l" . As k is a flat function of /?, is not crucial the 
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Lattice sizes 


Ul{Pc[uj = 0.5)) 


Ul[P,{lo = 1)) 


Ul{I3,{^ = 2)) 


6 


0.4435(15) 


0.4437(12) 


0.4438(12) 


8 


0.4406(15) 


0.4409(12) 


0.4413(12) 


10 


0.4407(14) 


0.4411(16) 


0.4414(14) 


12 


0.436(4) 


0.437(3) 


0.438(3) 


16 


0.435(3) 


0.436(3) 


0.437(3) 


20 


0.429(5) 


0.431(5) 


0.433(5) 


24 


0.428(6) 


0.430(7) 


0.433(7) 



Table 2.2: Binder cumulant for various lattices sizes at the extrapolated 
critical point for lo — 0.5,1,2. 



point where we actually measure. The results displayed in table 2.3 have 



been obtained measuring at the crossing point of the Binder parameters for 



lattice sizes L and 2L using ( 2.28 ) 



For measuring a/v we study the scaling of the specific heat 

C.^, ,2,37) 

We expect that C scales as A + BL°'^'^, where A is usua lly non-negligible. 
In Figure 6 we plot the specific heat measuring at ( ^.35 ), as well as at the 



peak of the specific heat, as a function of L. We observe a linear behavior 
for intermediate lattices. For the largest lattice the slope decreases. The 
weak first order behavior (a/i^ = 1 for small lattices that becomes d 
for large enough sizes) seem hardly compatible with our data. If we neglect 
the A term (what is asymptotically correct), and compute the exponent 
using eq. ( ^.28| ) we obtain a/v « 0.3 for intermediate lattices that reduces 
to a/iy = 0.15(2) for the (20,24) pair. 

However, to give a conclusive answer for the value of a statistics on 
larger lattices are mandatory. 

Magnetic Critical exponents: 7, /3 

The exponents 7 and P can be obtained respectively from the scaling of 
susceptibility and magnetization: 

X = V{M^) - r'/" (2.38) 

M - L-/^/" (2.39) 
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Figure 2.7: Specific heat at the peak (triangle symbols) and at /3 
(cross symbols) as a function of the lattice size. 



-0.7068 



Where M is an order parameter for the phase transition. In Figure 7 
upper part, we plot the quotient between Mpaf for lattices L and 2L as 
a function of the quotient between the Binder cumulants for both lattice 
sizes. 

For large L in the critical region we should obtain a single curve, the 
deviations corresponding to corrections to scaling. In the lower part of 
Figure 7 we plot the same function for susceptibility. 

The values for 7 and f3 are summarized in Table |2.3|. 



2.6.3 Logarithmic corrections 

We now address the question of the possibility of logarithmic corrections 
in the AF 0(4) model. For the thermal critical exponents, the situation 
seems clear, they are compatible with the classical exponent 1^ = 0.5. For 
the magnetic exponents, the situation is more involved. In principle, one 
can think that they disagree from MFT due to logarithmic corrections. We 
have no perturbative predictions about the form in which these corrections 
would affect for the AF case. However, one expects that such corrections 
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Figure 2.8: Quotients to obtain (3/v and 'y/v 



slightly modify the critical exponents, as occurs in the FM case. It could 
be possible that logarithmic corrections modify largely the previous critical 
exponents and drift them to the FM ones. To sort this out, we have consid- 
ered the possibility of a behavior FM like, so that ^ L(lnL)^/^. In the 
lower part of Table 2.3 we quote the values of the critical exponents for the 
PAF phase transition when logarithmic corrections are included (formula 
(2.32)). We see how in effect the magnetic critical exponents are too far 
from the classical ones for being the result of a logarithmic correction to the 
MFT predictions. It is interesting to compare this situation with that in 
the RP^ model in c? = 4 where small deviations from MFT exponents 
can be explained as logarithmic corrections. 



2.7 Conclusions and open problems 

We have studied the phase diagram of the four dimensional 0(4) model 
with first and second neighbors couplings. We have found that the pres- 
ence of couplings with opposite signs makes frustration to manifest in form 
of antiferromagnetism in different dimensionalities. Whether or not these 
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Lattice sizes 




aiv 


V 


(without logarithmic corrections) 


6/12 


2.417(3) 


0.791(4) 


0.474(10) 


8/16 


2.403(3) 


0.792(6) 


0.483(8) 


10/20 


2.410(2) 


0.790(4) 


0.471(6) 


12/24 


2.403(5) 


0.797(5) 


0.483(7) 


20/24 


2.398(5) 


0.802(4) 


0.487(6) 


(with logarithmic corrections) 


6/12 


2.301(3) 


0.849(4) 


0.484(9) 


8/16 


2.300(3) 


0.850(5) 


0.489(7) 


10/20 


2.315(2) 


0.843(3) 


0.488(5) 


12/24 


2.314(2) 


0.842(5) 


0.487(5) 


20/24 


2.317(5) 


0.839(4) 


0.498(5) 



Table 2.3: Critical Exponents for the PM-PAF phase transition in the F4 
lattice 



vacua have interest from the point of view of Quantum Field Theory is an 
open question. The transition between the paramagnetic phase, and the 
AF plane ordered one is compatible with very weak first order and with 
higher order. To shed some light on this question we have studied the evo- 
lution of the effective critical exponents in the limit /3i = 0. We found that 
up to L = 24 the exponents are in disagreement with the Mean Field pre- 
dictions. Specifically, from our 7/1/ estimation (or jv using hyper-scaling 
relation) the exponent 77 asociated with the anomalous dimension of the 
field is 77 « —0.4. This fact itself would imply the non-triviality of the 
theory because Green functions would not factorize anymore. We have also 
measured the Binder cumulant at the critical point, finding that it stays 
almost constant when increasing the lattice size. 

However, before claiming about the important issue of a non-trivial 
phase transition in four dimensions we have tried to acomodate the data in 
two possible scenarios: 1) logarithmic triviality; 2) weak first order scenario. 

The mere observation of critical exponents different from the Mean Field 
ones is in principle still compatible with a logarithmic trivial theory. We 
have contrasted our data with the perturbative predictions derived in the 
FM case, which do not fit the results. It would be pos sible to obtain 
triviality also with a logarithmic exponent in equation ( 2.32| ) different from 
1/4. We can fix the critical exponents to its mean field value and compute 
this parameter from the numerical data. The results obtained shown a non- 
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Figure 2.9: Effective j/i/ exponents in 0{N) witli N=2,3,4. 

asymptotic behavior with values ranging from 0.8 to 1.2 for the lattices 
used. In short, this model encounters the typical problems observed by 
other triviality studies, which make the numerical approach to this issue 
so difficult: one would need much stronger analytic insight to rule out 
logarithmic triviality. 

In the other hand, the possibility of being behind a very weak first 
order phase transition is a possibility that cannot be discarded a priori. 
Our data do not evidence any two state signature. Moreover, the stability 
of our measure of j/v for lattices ranging from L = 6 to L = 24, which 
are more than a hundred standard deviations apart the MF value, makes 
this hypothesis rather unlikely. However the question is wether or not this 
hypothesis can be ruled out beyond any doubt. The phase transition in 
this region for the Ising case, 0(1), has been shown to be first order [p^ . 
Lately, the models 0(2) and 0(3) have been considered with the purpose 
of checking wether or not the order of the phase transition changes with N 
in lattice sizes as larger as L = 48. Evidence is given for a first order 
character of the transition in 0(2). In 0(3) the transition is weaker, and 
despite no double peak structure is observed, the behavior of the effective 
exponents leaves little doubts about the first order character of 0(3) as well 
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(see Figure 2.9). In the 0(4) case we cannot observe any particular trend 
in the effective ^jv . Indeed, since the transition seems to be continuous 
in the hmit — > oo, = 4 might be a sort of critical N value from which 
on the transition is second order. 

However, from the observations in 0(2) and 0(3) we would be prone to 
consider the weak first order scenario as the most plausible. If the transition 
turns out to be first order, it would be weaker than the observed in 0(3), 
and hence lattices larger than L — 48 would be needed to observe such 
asymptotic behavior in the critical exponents. 
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Chapter 3 

La transicion de fase en el 
modelo SU(2)-Higgs 
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3.1 Introduccion 

Los modelos Higgs describen la autointeraccion del campo de Higgs y su 
acoplamiento con los campos gauge SU(2)i, y U(l)i'. 

En la aproximacion en la que se desprecian los acoplos de Yukawa entre 
el campo Higgs y los fermiones, estos modelos estan considerados como 
una razonable descripcion del sector electrodebil del Modelo Estandar. En 
particular se persigue mediante su formulacion en la red dar cuenta de los 
aspectos no perturbativos de la teoria electrodebil, como por ejemplo la 
generacion de masa mediante el mecanismo de ruptura espontanea de la 
simetria. 

Los acoplamientos gauge dc U(l)y y SU(2)i se denotan respectivamente 
por g'/2 y g. La relacion entrc olios as a traves del angulo de Weinberg, 
concretamcntc la relacion es g'/g = tan^w, puesto que este angulo es 
pequefio (sin^ 6'w ~ 0.23), podemos despreciar en un primer paso el grado 
de libertad U(l)y. En este sentido SU(2)-Higgs puede considerarse como 
el modelo Higgs mas simple. 

El interes del estudio de estos modelos no es tanto la necesidad de 
determinar parametros medibles en la iiaturaleza, sino mas bicn profundizar 
en el entendimiento teorico de la formulacion de la teoria electrodebil. A 
diferencia de la QCD, esta teoria no es asintoticamente libre, el punto fijo 
gaussiano a acoplamic;nto cero es inestable ultravioleta, y la existencia de 
un punto fijo no trivial se hace necesaria para la formulacion de la teorfa 
en el continuo. 

La accion del modelo SU(2)-Higgs consiste en una parte pura gauge mas 
un termino de interaccion escalar-gauge: 



La interaccion entre la parte escalar y la gauge esta controlada por el 

hopping parameter, k. Una normalizacion adccuada cuando se esta intcrc- 
sado en trabajar en el limite de autoacoplamiento cuartico infinito (regimen 
no perturbativo) es la siguiente: 



con esta normalizacion el Hmite A ^ oo equivale a fijar el modulo del 
campo Higgs a la unidad. 

La parte pura gauge, [U] , se expresa como la suma sobre las plaquetas 
orientadas positivamente: 



S[U,^]x = Ss[U] + S^[U,^] . 



(3.1) 




(3.2) 
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5g=/3^(l-itrC/p) , 
p 



(3.3) 



donde (3 = A/g^. En el Ifmite /3 ^ oo el campo gauge se desacopla y 
se recupera el modelo el cual presenta una transicion de fase entre una 
region simetrica (pequeno n) y una region en la que la simetri'a se ha roto 
espontaneamente (k grande) 

Cuando /? es finita (g ^ 0) esta transicion de fase se prolonga en una 
Imea de transiciones de fase en el interior del espacio de acoplos HJ ^, D, 
como se muestra en la figura |3.l| . 

La continuacion de la fase simetrica de cuando el acoplo gauge es 
distinto de cero es una region conGnante en el sentido de que los campos 
gauge SU(2) confinan las parti'culas escalares de materia. En este sentido 
la fase simetrica de SU(2)-Higgs es analoga a QCD. 

La continuacion de la fase de simetrfa rota cuando el coupling gauge 
es distinto de cero es la fase Higgs. En esta fase, mediante el mecanismo 
de ruptura espontanea de la simetrfa, el campo de Higgs adquiere un valor 
esperado en el vacfo distinto de cero. Los bosones vectoriales adquieren su 
masa a traves de su acoplamiento con el campo escalar. Los tres bosones 
de Goldstone del modelo se hacen masivos al mezclarse con los bosones 
gauge vectoriales. El resultado es un isovector con spin 1, el boson W. 

Esta Ifnea de transiciones de fase termina en un valor finito del espacio 
de acoplamientos, estando las fases Higgs y confinante conectadas analiti- 
camente. 

En cuanto al orden de la transicion de fase, para pequenos acoplamientos 
existen argumentos perturbativos debidos a Coleman y Weinberg |^ apun- 
tando que la transicion deberi'a ser de primer orden. El estudio numerico 
de la transicion de fase para pequeiios e intermedios valores de A parece 
confirmar esta conjetura H, ||]. 

Sin embargo se observa que la transicion se debilita (las discontinuidades 
son mas pequefias) conforme A aumenta. En el Ifmite A = oo la opinion 
generalizada es que la transicion sigue siendo de primer orden aunque ex- 
tremadamente debil, si bien no existen evidencias que justifiquen esta creen- 
cia. Los liltimos resultados numericos fueron obtenidos por Montvay en 
1985, en redes no demasiado grandes y a pesar de observar metaestabili- 
dades en los ciclos de histeresis no se pudieron obtener resultados con- 
cluyentes [ ]Tl[ . 

Tratando de arrojar algo de luz sobre este punto hemos enfocado el 
problema desde un punto de vista distinto, formulando el modelo SU(2)- 
Higgs en un espacio de parametros extendido. Sin embargo la motivacion 
no es solo hacer un estudio con alta estadfstica del modelo, sino tratar de 
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Figure 3.1: Grafico esquematico del diagrama de fases del modelo SU(2)-Higgs 
con modulo fijo (A = oo). 

extraer las propiedades generales que presentan las transiciones de primer 
orden debiles en la red en modelos gauge acoplados a campos escalares. 
Con este proposito se ha anadido a la accion estandard un acoplo extra 
entre el campo gauge y el campo escalar conectando segundos vecinos en 
la red. Este acoplamiento nuevo sera usado como un parametro para re- 
forzar o debilitar la transicion, permitiendo asi estudiar su mecanismo de 
debilitamiento cuando este acoplamiento se hace tender a cero. 
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Abstract 



The properties of the ConRnement-Higgs phase transition in the SU(2)- 
Higgs model with fixed modulus are investigated. We show that the system 
exhibits a transient behavior up to L—24 along which, the order of the phase 
transition cannot be discerned. To get stronger conclusions about this 
point, without going to prohibitive large lattice sizes, we have introduced 
a second (next-to-nearest neighbors) gauge-Higgs coupling (K2). On this 
extended parameter space we End a line of phase transitions which become 
increasely weaker as K2 — > 0. The results point to a first order character 
for the transition with the standard action (k2 — 0)- 

3.2 The Model 

The SU(2) lattice gauge model coupled to an scalar field, in the fundamental 
representation of the gauge field can be described by the action 

Sx = /3 ^[1 _ i tr C/p] - ^K, *^(a;);7^(x)*(x + fi) + 

+A^[*^(a;)*(a;) + ■ (3-4) 

X X 

Where U^(a;) represents the link variables, and Up are their products 
along all the positive oriented plaquettes of a four-dimensional lattice of 
side L. 

The scalar field at the site x is denoted by being A the parameter 

controlling its radial mode. In the limit {X — 00, f3 — 00) the model 
becomes a pure 0(4)-symmetric scalar model. 

As we pointed out in the previous section, the PT line ends at some 
finite value of the parameters (/3, ki) The endpoint moves towards larger /3 
values as A increases. For A > 0.1 the endpoint crosses to the /? > region, 
and in the limit A = 00 the phase transition ends at (ki ~ 0.6, f3 ~ 1.6) 
It is commonly believed that the transition at this point is second order with 
classical critical exponents , however a careful numerical study would be 
necessary. 

In the scaling region, and for finite A, the phase transition turns out to 
be first order. Also, the transition becomes weaker as /3 or A increases. In 
particular in the limit (3 = 00 (spin model) the transition is second order 
with classical critical exponents ||l|. 

In the limit A = cxd the situation is less transparent, and it is not clear 
whether the transition is weak first order or higher order. 
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The study of the model with the parameter A = cx) is equivalent to fixing 
the modulus of the Higgs field, • 4> = 1. In the pioneer work, [|oj , the PT 
was considered first order for finite A and second order for A = oo. Later, 
larger statistics, and the hypothesis of a universal behavior of the PT for all 
values of A, seem to point to a first order character of the PT in the scaling 
region. Nowadays, though it is generally believed that the transition is still 
first order in this limit, the numerical proofs |^ on which rely these 
statements are not conclusive, as the authors safely conclude, because the 
statistic and lattice sizes are not enough for excluding the possibility of a 
higher order phase transition in the edge X — oo. 

The model with fixed modulus is described by the action 

Soo=PY.^l-^ti-Up]-^KiY,^i-^\x)U^{x)^{x + pi) . (3.5) 

As we shall show below, we have simulated this model up to lattices 
L—24, and the result is still compatible with a second or higher order PT, 
however, we have no indications of asymptoticity in the behavior of the 
observables, and the results are compatible with a very weak first order PT 
too. This means that larger lattices are needed to overcome these transient 
effects, but the added difficulty here is that such lattices would suffer of ter- 
malization problems, and severe autocorrelation times. Altogether makes 
this approach too CPU expensive for nowadays computers. 

In order to get a more conclusive answer without going to prohibitive 
large lattice sizes, we have studied the model in an extended parameter 
space. For this purpose we have introduced a second coupling between 
the gauge and the scalar field connecting next-to-nearest neighbors on the 
lattice, in such a way that the new action reads 

(3.6) 

Within this parameter space we expect to get a global vision on what the 
properties of the PT are, and also, to give a stronger conclusion about the 
order. In the region of K2 positive, (competing interaction effects appear if 
K,2 < 0) this extended model is expected to belong to the same universality 
class than the standard one (^2 = 0) since both models posses the same 
symmetries. The effect of the new coupling K2 > is to reinforce the 
transition but should not change the order if the system has not tricritical 
points. An example of such behavior appears in the 0(4)-symmetric a 
model with second neighbors coupling. This model presents a {k^, Kj) line 
of phase transitions which is second order, since the model with K2 = 
shows a second order PT too 
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The action (3.6) has the foUowing symmetries: 



• Kl = K2 = 0. 



(3.7) 



Ki — > — Ki, ^{x) fixed 

Ki —Ki, U^{x) fixed 

*(a;) ^ (-1)^^=0'''' *(a;) 



(3.8) 



(3.9) 



The action is not symmetric under the change K2 — > — K2. The exis- 
tence of couphngs Ki and K2 with opposite signs would make frustration 
to appear, and very different vacua are possible In this work we are 
interested in the regions free of frustration effects. Taking into account the 
symmetry properties of the action, the phase diagram in the region ^2 > 
will be symmetric with respect to the axis ki = 0, and hence, we can 
restrict the study to the quadrant ki > 0. 

We define the normalized energy associated to the plaquette term 

i?o = ^^(l-itrC/p) , (3.10) 

p 

and also the energies associated to the hnks 

£^i = ^E^t^*^(^)^M(2^)*(2^ + M) > (3-11) 

E2 = ^J2 l^r^Hx)[U^{x)U,{x + ti) + U,ix)U^ix + ,^)Mx + fi + iy) 

(3.12) 

where Ni„ = 6V, Ni^ = 12V and Ni, = 4V. 

With these definitions Eq when P 00 and iJ^ — s- 1 when — > cx). 

On the three-dimensional (/3, ki, K2) parameter space we consider the 
plane /3 — 2.3. On this plane there is a PT line (k^,K2). We expect to 
learn on the properties of this PT on the region (k^ 7^ 0, K2 ^ O), where the 
signals are clearer, with regards to applying what we learn to the standard 
case K2 = 0. 
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Figure 3.2: Phase diagram obtained from the MC simulation. 



To monitorize the strength of the phase transition, we measure the 
existence of latent heat, and the behavior of the specific heat. As we shall 
see, for ki = the transition is first order, with a clearly measurable 
latent heat. We will see how this transition weakens along the PT line for 
increasing values of ki. 



3.3 Numerical study 

We have simulated the model in a lattice with periodic boundary con- 
ditions. For the update we have employed a combination of heat-bath and 
over-relaxation algorithms (ten over-relax sweeps followed by a heat-bath 
sweep). For the simulation we used the RTNN machine, consisting of a net- 
work of 32 PentiumPro 200MHz. processor. The total CPU time employed 
has been the equivalent of 3 years of PentiumPro. 

Monte Carlo methods provide information about the thermodynamic 
quantities at a particular value of the couplings. We have used the Spectral 
Density Method (SDM) to extract information on the values of the 
observables in a finite region around the simulation point. In particular it is 
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0.380 0.385 0.390 0.395 

Figure 3.3: Normalized energy distribution for L=12, A = oo, ki=0.395, 13—2.3 
from (dotted line), compared with the distribution we obtain at the same 
couplings (solid line) in L=12 too, when statistics is increased by one order of 
magnitude. 

useful to have a precise location of the coupling where some observables have 
a maximum, as well as an accurate measure of the value of that maximum. 

From the Monte Carlo simulation at some coupling k, we got the his- 
togram H{E) which is an approximation to the density of states. Using the 
SDM approximation the probability of finding the system with an energy 
E at a different coupling k! can be written as: 

P^\E) cx H{E)e'^'''-''^^'^ . (3.13) 

The region of validity of the SDM approximation is Ak ^ 1/{V(j), being 
cr the width of the distribution H{E). Although a gets the maximum values 
in the critical region, the approximation has been very useful, specially for 
tuning the couplings where to measure. 

Concerning the lattice sizes, we have used lattices ranging from i=6 
to L=24. For the small lattices {L—6, 8 and 12) we have done 4xl0^r 
iterations, being t the largest autocorrelation time for the energy, which 
ranges from r « 10 in L=8 to r « 35 in L=24. For the largest lattices. 
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L=2Q and 24, we run up to lO^r MC iterations. 

The statistical errors are computed with the jackknife method. 

3.4 Results 

We shaU make the discussion with the first-neighbors hnk energy, Ei, but 
as far as the critical behavior is concerned, we could carry out the analysis 
with any of the energies. We remark that an appropriate linear combination 
of El, E2 and Eq could give slightly more accurate results [Q. 

We have considered fixed values of ki (0, 0.02, 0.1, 0.2 and 0.3) and 
sought the K2 critical for every line, K2(ki). We have also studied the case 
K2= varying ki which corresponds to the usual SU(2)-Higgs model. 

The SDM has been used to locate the apparent critical point, defined 
through the specific heat behavior. From the specific heat matrix: 

Cl'^L)^^, (3.14) 

we obtain the best signal for C^'^(i), which can be calculated as (we shall 
omit the superscript from now on) 

a(i) = 4L''{{EiE2) - {Ei){E2)) . (3.15) 

In a first order phase transition, G^^^^{L) behaves, asymptotically, pro- 
portional to the volume, L'^. If the PT is second order, the dominant 
behavior for Cy{L)"^'^'^ is L"/"^ which diverges too provided that a > 0. At 
the upper critical dimension a = 0, and one has to go further the leading 
order, appearing logarithmic divergences p^ . 

As a consequence of this divergent behavior, in a finite lattice Ci,(L) 
shows a peak at some value of the coupling which will be taken as apparent 
critical point, n'^ jL). 

In Figure |3.2| we plot the critical line (k^, K2)- This line is obtained by 
extrapolation to the thermodynamic limit according to ^2(00) = K2{-^) ~ 
AL^'''. We point out that this extrapolation is valid only in first order 
phase transitions. If the PT is second order the power in L is {—1/v). 

In continuous PT scale invariance holds, and the thermodynamic mag- 
nitudes, such as the specific heat, or susceptibilities do scale. However if 
the transition is first order the correlation length remains finite and hence 
there is not scaling properties, and no critical exponents can be defined. 

Nevertheless in a first order PT we can ask how large is the lattice 
we need in order to observe the asymptotic behavior of C^. In particular, 
with abuse of language one can measure a pseudo a/v exponent to get 
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Figure 3.4: Ei distribution on the axis ki = on a L=8 lattice at {k2 — 0.15, [3 — 
2.3). The cosine fit is very accurate in spite of the finite (3 value, in this region of 
parameters the pure gauge term couples slightly both sub-lattices. 

insight on the nature of the PT: the larger is the lattice we need to measure 
= the weaker the PT is. Following this, first order PT can be 

classified according to their degree of weakness. 

However, the so called weak first-order PT appear often in literature (see 
[ p^ , and references there in) as PT characterized by a transient behavior 
with a non-measurable latent heat. Let be the correlation length of the 
system at the critical point in the thermodynamic limit. In a finite lattice of 
size L, the first order behavior will be evidenced if i > For lattice sizes 
much smaller than the system will behave like in a second-order PT, since 
the correlation length is effectively infinite. As an example, in Figure |3.3| 
we plot the energy distribution in a L=12 lattice at {f3 = 2.3, m — 0.395) 
obtained in pT| ], compared with the one we obtain at the same couplings 
and in the same L, when the statistics increases by one order of magnitude. 
We observe that the order of the PT can not be discerned at this volume, 
even when the statistics is enough. Termalization effects can also contribute 
to mistake the histogram structure. 

The entire line {ki,K2) is first order, but the weak character increases 
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as K2 — > 0. We will make a quantitative description of the weakening 
phenomenon by studying the specific heat, and the latent heat. 

But before going on, we shall make a remark concerning the behavior 
on Ki = 0. As we pointed out, the system is symmetric under the change 



Ki — !■ — Ki. The transformation (3.9) maps the positive ki semi-plane with 
energy Ei, onto the negative ki semi-plane with energy -Ei. The transition 
across this axis is first order because the energy is discontinuous. In the 
limit /3 ^ CO, and in ki ~ 0, the system decouples in two independent 
sublattices, each one constituted by the first neighbors of the other. The 
first neighbors energy for this system is proportional to cos 9, being 6 the 
angle between the symmetry breaking direction of the scalar field in both 



sub-lattices. In Figure 3.4 we plot the Ei distribution for L=8 at (ki = 
0,K2 = 0.15) and /? = 2.3. We see that the agreement with a cosine 
distribution is quite good in spite of the finite /3 value. 
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Figure 3.6: MC evolution of Ei at (ki = 0.2, n^L)) for L=8 (lower part), L=12 
(middle) and L—16 (upper part). 



3.4.1 Latent Heat 



Along the apparent critical line wc have done simulations for different lattice 
sizes and stored the plaquette and links energies to construct the histograms 
for the energy distributions. In a first-order phase transition the energy has 
a discontinuity which manifest in the appearance of latent heat, AE. This 
quantity is not well defined in a finite lattice, so we measure the distance 
between the two maximum of the energy distribution, and extrapolate to 
the thermodynamic limit. The drawback of this approximation is that the 
maxima of the energy distribution are difficult to discern, since this function 
at the apparent critical point is very noisy We have used a cubic spline at 
the maxima in order to get a more reliable estimation (Figure 5). 

In Figure we show the MC evolution of Ei for L=8, 12 and 16 at 
Ki = 0.2. In L=8 the latent heat is not clearly measurable. We observe in 
the MC evolution how th e tw o-state signal becomes cleaner as the lattice 
size increases (see Figure 3.7). 

In Figure 3.7 we plot the distribution of Ei at ki = 0.2 at the apparent 
critical point k,2{L) for L=6, 8, 12 and 16. A remarkable stability of AEi 
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Figure 3.7: Normalized distribution of Ei at (rei = 0.2, K2(I/)) for L=6, 8, 12 
and 16. 



with the volume is observed. This is a common feature for all values of ki. 

As we have already pointed out, the transition weakens when increasing 
Ki and larger lattices are needed in order to observe a measurable latent 
heat. We give a quantitative description of this fact in Figure 3.S, where 
the distribution of Ei for several values of (ki,K2(12)) is displayed. The 
two-state signal is no longer measurable in L=12 at ki = 0.3. The first 
evidences of two-states appear in L—20 (see Figure |3.9| ) . but from its energy 
distribution we can only give an approximate value for AEi (L — 20) since 
the two-peaks appear too close to each other. 

In Figure 3.10 we plot AEi{L) and AE2{L) as a function of l/L*, in 
order to get AEi in the thermodynamic limit with a linear fit. Finally we 
quote these values in Tabic |3.l| , together with the change in the action ( ^ ) 
between the two phases. 

From the energy distributions at K2 = 0, see Figure 3.11 we have no 
direct evidences of the existence of latent heat. However, on the larger 
lattices one can observe non-gaussianities in the energy distributions. Such 
asymmetries could precede the onset of clear two-peak structures in larger 
lattices, however this is just a guess. We conclude that no information con- 
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Figure 3.8: Normalized distributions of Ei for rei = 0.02, 0.1, 0.2 and 0.3 at 
^2(12). 



cerning the order of the PT can be obtained from the energy distributions 
up to L=24. 



3.4.2 Specific Heat 

We have done MC simulations at the points predicted by SDM, K2{L)), 
in order to measure accurately the peak of Cv{L). 

As an example we show in Figure 3.12| the value of Cy (L) around its 



maximum for various lattice sizes at ki — 0.2 (upper plane) and at ki 
0.3 (lower plane). We observe that the maximum of the specific heat grows 
slower when increasing ki, indicating a weakening in the PT. 

In Figure |3.13 we draw C™'*'^(L) relative to C^'^^{6) as a function of 



the lattice size. The values have been normalized to C™'*'^(8)/C™'*'*^(6) in 
order to compare distinctly the behaviors for different ki . The slope of the 
segment joining the values of C^^^^^ in consecutive lattices gives the pseudo 
ajv exponent. We observe that such slope is approximately 4 at ki = 
0.02, 0.1 and 0.2 for all the volumes we compare. However the transition 
at K\ = 0.3 evidences much more weakness. We do not have evidences 
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Figure 3.9: Normalized distribution of Ei at Ki = 0.3 for L=12,16 and 20 



of asymptoticity in G™'"^^ till L=20, as could be expected from the energy 
distributions (Figure [3^ ). The slope of the segment joining C™'*''(16) with 
Qmax^20) is 3.05(12) which is almost the asymptotic value expected for a 
first order PT. 

As expected, only if the two-peak structure is observed in the energy 
distributions and l^E is stable, the maximum of the specific heat will grow 
up like the volume, L'^. 

At K2 = we are within the transient region even for L=24. We remark 
that in this case is defined by the element C^'"'^ of the specific heat 



matrix (3.14) since this is the most natural choice at this point, and also 



is the best signal we measure. As we observe in Figure 3.13, at K2 = 
(-imax^^-j seems to tend to a constant value as oo up to L—20. From our 
previous discussions we should conclude that either the correlation length 
at the transition point is much larger than the lattice size up to L=20, 
or the transition is second order with a = in the thermodynamic limit. 
However, in L = 24 things are changing, C^'^^{L = 24) starts to run away 
of this quasi-plateau, and the pseudo ajv exponent grows again. As can 



be observed in Figure 3.13, at k.\ = 0.3 the pseudo ajv exponent decreases 



for the segment L = 12-16, with respect to the value in the segment L = 



11 



coupling 




Ai;2(oo) 


AS' 


Kl = 




0.0366(8) 


0.0134(12) 


Ki = 0.02 


0.0162(6) 


0.0347(5) 


0.0137(13) 


Ki = 0.1 


0.0162(7) 


0.0345(9) 


0.0094(10) 


Ki = 0.2 


0.0179(7) 


0.0201(8) 


0.0078(12) 


Ki = 0.3 


RiO.006 


w 0.012 


w 0.0026 



Table 3.1: A£(oo) for Ei and E2, and variation of the action. 
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Figure 3.10: ABi (upper plane) and Ai52 (lower plane) as a function of 1/L*. 
The two-peak structure is not clearly observed in L=6 at any ki value. The 
values quoted for this lattice size are upper bounds. 



8-12. The lattice L = 20 is enough to overcome the transient region, but 
the behavior is qualitatively the same that in K2 = 0, though the transition 
is stronger. 

We believe that this behavior is general for weak first order transitions 
in four dimensions. There exists a transient region in which the correlation 
length is effectively infinite compared with the lattice size, and the system 
behaves like suffering a second order PT with thermal index a w in the 
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Figure 3.11: Normalized distribution of Ei at K2=0, in L=16, 20 and 24 at the 
peak of the specific heat 



thermodynamic limit. 

3.4.3 Binder Cumulant 

In order to check the consistency of our results we have also considered the 
behavior of the Binder cumulant 

V.^l-^, (3,16) 

This quantity behaves differently depending on the order of the PT. If 
the transition is second order the minimum of the cumulant, V^™^ approach 
2/3 in the thermodynamic limit. However if the transition is first order, 
V^"'" tends a value smaller than 2/3 reflecting the non-gaussianity of the 
energy distribution at the transition point. 



In Figure 3.14 we plot V™™ for several ki values. 

For those values of ki in which the PT is distinctly first order, the 
minimum of the Binder cumulant stays safely away from 2/3, as we observe 
y^nin extrapolated to oo is 0.65401(4) at k\ = 0.1. However, this value 
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Figure 3.12: C™^''(L) at m = 0.2 (lower plane), at m = 0.3 (middle) and at 
K2 — (upper plane). The dotted line is the SDM extrapolation. 



reachs 0.66637(5) at ki = 0.3, and 0.66657(8) at K2 = 0. Again we find a 
tight difference between a very weak first order PT and a continuous one. 

For the safce of discussing quantitatively the order of magnitude of the 
latent heat in the hmit K2 = 0, from the energy distributions we find that 
at K2 = one can approximately locate one the peaks of the energy at « 
0.273. The other should be at certain Ef, = Eq — A, being A the latent 
heat. 

In the thermodynamic limit the energy distributions are two delta func- 
tions situated at E^ and E^, then 



rain 
oo 



HEl 



El? 



(3.17) 



If we use V™" 



— 0.66657 and Ea in ( 3.17 ) the value we got for the latent 
heat is A « 0.006 which is of the same order as the one expected from the 
histograms. 
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Figure 3.13: C™^ for the various ki values and K2 = 0. We have normaUzed the 
values with respect to C™'^(8)/C™'"(6). The slope of the segments is indicated 
when smaller than 4. 



3.5 Conclusions 

The order of the Confincmcnt-Higgs phase transition in the SU(2)-Higgs 
model with fixed modulus is a highly non trivial issue. We have used an 
extended parameter space, in order to get a global vision on the problem. 
On this extended parameter space we have found a line of first order phase 
transitions which get weaker as k;2 ^ 0. We have also observed that, 
because of the computer resources needed, it is too ambitious trying to 
measiirc two-peak energy distributions in the limit K2 = 0. However, on 
this point, we can extract conclusions from the behavior of the specific heat. 

As we have discussed along the paper, a fake second order PT seems 
to exists for a range of L in very weak first order phase transitions. We 
have applied Finite Size Scaling properties along this transient region to 
compute a pseudo ajv critical index. We want to be extremely careful at 
this point, this computation is completely meaningless when the transition 
is first order, since Scaling does not hold, but it can be used as a technical 
tool to catalogue the PT when there is no direct evidences, as in this case. 
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Figure 3.14: VT'" for m = 0.1, 0.2,0.3 and k2 = 0, as a function of 1/L''. The 
dashed line represents the value 2/3. 



Using the relation a — 2 — i/d, we got i/ varying in the interval (0.36, 0.41) 
in the range L = 8, . . . , 20. Calculated from i=20 and i=24, « 0.35. We 
expect this behavior to be transitory, and when going to larger lattices sizes, 
if the transition is second order, v should reach its mean field value v = 1/2. 
If the transition is first order this value should go to 1 /d, indicating that the 
specific heat maximum grows like the volume L"^. We believe that this is 
the case, since the pseudo v exponent in i=:24, instead of approaching 1/2, 
starts to decrease. An example of weak first order PT showing a similar 
behavior is described in [igt] . 

In what concerning the motivation of introducing a second coupling, we 
pointed out that K2 should not change the order of the PT because do not 
change the symmetry properties. This argument is heuristic, but the phase 
diagram we found supports this assertion. As far as the order of the PT is 
concerned, we think that this approach can be useful when dealing with PT 
of questionable order in the sense that it is not clear whether the transition 
is weakly first order or higher order. The hope is that it could be applied 
to other more controversial models. 
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Chapter 4 



El modelo U(l) gauge en 
D=4 con topologias 
toroidal y esferica 
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4.1 Introduccion 



La transicion de la teoria libre de Dirac a la electrodinamica cuantica (QED) 
implica reemplazar una simctria global por una local es decir, hacer gauge 
la teoria. Concretamente de la accion correspondiente a la teon'a libre de 
Dirac: 

4°) = J dx^^{x){ird^ - Af)V(x) , (4.1) 

invariante bajo transformaciones globales de U(l), es decir del tipo G = e^^, 
se pasa a la electrodinamica cuantica (QED) haciendo que la invarianza 
global sea tambien local, esto es, invariante bajo cambios locales de fase 
pertenecientes a U(l), es decir G{x) = e^^^^\ Esta transformacion implica 
introducir un potencial gauge Af^{x) y reemplazar las derivadas ordinarias 
por derivadas covariantes: — + ieA^{x). 

Sin embargo para que el campo gauge juegue un papel en la dinamica 
hay que afiadir a la accion original un termino cinetico que permita su 
propagacion. Para ello se introduce el tensor T^^'- 

^M-(x) = -'-[D^,D,] = d^x) - duA^[x) , (4.2) 

de manera que la accion que describe la QED en el espacio euclideo 
queda: 

ccont c I C 

'^QED - + JF 



dx^Tf,,{x)Tf,,{x)+ / dx^i}{x){l-1^'D^, + M)^{x){A.2,) 



1 



En lo que sigue estaremos interesados en el estudio de la parte pura 
gauge en la red, donde una posible formulacion de la teoria del continuo 
esta dada por la accion de Wilson §: 

donde U^v{n) es el usual producto ordenado de las links que forman 
cada plaqueta 

[/^.(n) - U^,{n)U,{n + ti)Ul{n + v)Ul{n) . (4.5) 

En el limite en el que el espaciado de la red se hace tender a cero se 
recupera la accion del continuo, Sg, con la siguiente correspondencia: 
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U^Xn) = e*«'^^M(») ^ (4.6) 
= ^. (4.7) 
En efecto, usando la version discretizada de T^u'- 

^M-(n) = ^[{A,{n + A,{n)) - (A,,(n + u) - A^{n))] , (4.8) 
se tiene que la variable plaqueta se puede escribir como: 

U^^^in) « e'^'^'-^-(") , (4.9) 
y por tanto la accion en el li'mitc a —t queda: 

n 

El modelo U(l) compacto en d = 4 posee dos fases. Una es la region 
de acoplamiento fuerte en la cual los fotones estan confinados, y la otra 
es la region de acoplamiento debil en la que los fotones son portadores 
de un potential tipo culombiano. La transicion de fase que separa ambas 
regiones tiene lugar en /3 ~ 1 y ha sido ampliamente estudiada debido a 
la relevancia del modelo U(l) para estudiar QED. Sin embargo el estudio 
de esta transicion de fase ha resultado ser mas complicado de lo que se 
podia esperar en un principio, de su al menos aparente simplicidad. En la 
actualidad, las cuestiones relacionadas con el orden de la transicion de fase, 
y con los mecanismos que la provocan, siguen siendo temas controvertidos. 

La mayorfa de los estudios numericos han sido realizados en redes ciibicas 
en AD con condiciones de contorno periodicas, es decir en un hipertoro 
(TiT). Los estudios iniciales de este modelo apuntaban a que la tran- 
sicion era continua, puesto que las simulaciones numericas no evidenciaron 
ninguna serial de metaestabilidad [|[ §1 ■ Las redes mas grandes utilizadas 
en aquellos trabajos fueron L = 4, 5. Sin embargo, cuando los recursos de 
los ordenadores permitieron simular en redes mas grandes, se empezaron a 
observar metaestabilidades e histogramas con doble pico de L = 6 en ade- 
lante, revelando que la transicion tenia calor latente y por lo tanto serfa de 
primer orden |^, ||]. Esta idea fue corroborada por trabajos que aproxi- 
maron el problema desde el punto de vista del Grupo de Renormalizacion 

Actualmente, quedan un cierto numero de problemas abiertos. 
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1. Las simulaciones desde L — 6 hasta L — 12 revelan calor latente, pero 
su valor decrece con L. 

2. El cxponente cn'tico efectivo esta en el intervalo (0.29,0.32) que es 
distinto de lo que se espera en primer orden, v — 0.25, y tambien del 
valor correspondiente al segundo orden trivial, = 0.5. 

Estos hechos llevaron a considerar la posibilidad de que el calor latente 
fuese a cero en el Hmite termodinamico. Es decir, la transicion de fase en 
U(l) podri'a a pesar de todo ser continua si los dos picos que aparecen en 
las simulaciones fuesen un efecto de tamaiio finito [|l3j . 

Aparte de la accion de Wilson, se consideraron acciones extendidas para 
tratar de clarificar el problema mediante su formulacion en un espacio de 
fases extendido. En este contexto, una accion que incluye un termino 
proporcional a las plaquetas al cuadrado fue propuesta por G. Bhanot a 
comienzos de los '80 [l^ : 

5£w = /3EE[i-^(^A..+^M.)]+7EE[i-^(^'.+^^')] ■ (4.11) 

n /^<i^ n fj,<u 

Esta accion posee el mismo limite continuo que la accion de Wilson con 
la identificacion /3 + 47 = l/g^ . 

En este trabajo se estudio el espacio de parametros extendido (/?, 7) 
en L = 4, 5. Se encontro que la linea transicion de fase desconfinante se 
extiende para todo 7. En estas redes solo se observaron metaestabilidades 
en 7 > +0.2, sugiriendo que la linea de transiciones de fase Pdl) se hace 
de primer orden para 7 > +0.2, mientras que en 7 < +0.2 la transicion 
aparecia como de segundo orden. 

Simulaciones posteriores, siempre en redes TIT, ban mostrado la exis- 
tencia de calor latente para todo valor de 7, tanto positive como negativo, 
para reti'culos suficientemente grandes. Por tanto, el trabajo de Bhanot no 
es exacto cuantitativamente, pero cualitativamente sus conclusiones si son 
correctas: los valores de 7 positives refuerzan la transicion en el sentido que 
el calor latente crece cuanto mas positive es este parametro, mientras que 
los valores negatives de 7 debilitan la transicion y calor latente decrece. 

Actualmente se acepta con generalidad que la transicion es de primer 
orden para valores positives de 7. Sin embargo para 7 < sen posibles dos 
escenarios: 1) la linea de transiciones (3c{"f) es de primer orden para todo 7 
finito, positive o negativo; 2) Existe un 7tcp < punto tricrftico a partir 
del cual la transicion es de segundo orden |l3|| . 

Este problema fue estudiado desde el punto de vista del Grupo de Renor- 
malizacion [Q . Los resultados dan soporte al primer escenario descrito, en 
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concrete se encuentra un debilitamiento de la transicion al hacer 7 mas ne- 
gative, pero ello sin cambiar el caracter de primer orden, sino simplemente 
debilitandola. 

Hasta aquf se han descrito los resultados existentes en la topologi'a TiT . 

Sin embargo a lo largo del estudio de esta transicion de fase se han 
usado otras topologfas aparte de la toroidal. El motivo para ello es la 
idea debatida hace algunos afios sobre la influencia de los monopolos en 
la transicion de fase de U(l). Concretamente se argumento que debido a 
las condiciones de contorno periodicas existen monopolos no triviales, en 
el sentido de que no son contrafbles a un punto Al ser los monopolos 
configuraciones asociadas a mi'nimos locales de la accion, su existencia hace 
crecer la energia, y se sugirio que el salto en la energi'a observado en la 
transicion de fase se debe a la desaparicion de estos monopolos, que en la 
region desconfinada no existen. 

Estas hipotesis indujeron a trabajar en redes con condiciones de con- 
torno esfericas, puesto que en la esfera todos los monopolos son contraibles 
a un punto jl^, En concrete ha sido sugerido que los dos estados 

metaestables desaparecen cuando se trabaja en una red construfda con- 
siderando la superficie de un cubo en bD p^ , |l8| , que tiene la misma 
topologfa que la esfera. La ausencia de doble pico, junto la estimacion 
del exponente cri'tico u ~ 0.37 medido en 7 = 0, —0.2, —0.5 ha inducido a 
los autores en |Q a decir que los dos picos observados en el toro son un 
efecto de tamario finito, y que trabajando en la esfera se ve que la linea 
/3c (7) es critica para 7 < y esta caracterizada por un exponente critico 
1/3. 

El papel de los monopolos no triviales mencionados anteriormente parece 
estar descartado en cuanto a su influencia en el orden de la transicion de 
fase Por lo tanto a primera vista la situacion es bastante extrana. 

Puesto que que los monopolos no triviales no parecen jugar ningiin papel, 
uno esperarfa obtener antes el comportamiento asintotico en una red ho- 
mogenea e invariante por traslaciones, es decir la toroidal, que en una red 
manifiestamente no homogenea, como la superficie de un cubo en SZ^, y en 
la cual la invariancia traslacional se recupera solo en el Ifmite de volumen 
infinito. 

Como se ha venido exponiendo a lo largo de esta memoria, la ausencia 
de doble pico no prueba que una transicion sea de segundo orden. Ademas 
la aparicion de un exponente v ^ 0.37, intermedio entre el valor de primer 
orden y el de segundo orden trivial, ha sido encontrado en otros modelos 
como caracterizador del regimen transitorio de las transiciones de primer 
orden debiles @, l^]. 

Para tratar de arrojar algo de luz sobre este problema se ha hecho un 
estudio sistematico de la linea de transiciones de fase /3c (7) en el toro y en 
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la red con topologi'a esferica empleada en |16, Qq] y que denotaremos por 



JiS. Para esta red hemos empleado tamanos mas grandes de los medidos 
anteriormente para ver si la senal de doble pico aparece. En la red TiT 
hemos mejorado la estadistica que habia hasta la fecha, midiendo en redes 
hasta L — 2A para algunos valores de 7. 
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Abstract 



We have performed a systematic study of the phase transition in the 
pure compact U(l ) lattice gauge theory in the extended couphng parameter 
space (P, 7j on toroidal and spherical lattices. The observation of a non-zero 
latent heat in both topologies for all investigated 7 G [+0.2, —0.4], together 
with an exponent Vctf ^ when large enough lattices are considered, 
lead us to conclude that the phase transition is first order. For negative 7, 
our results point to an increasingly weak first order transition as 7 is made 
more negative. 

4.2 Description of the model and observables 

We shall consider the parameter space described by the extended Wilson 
action, which can be expressed in terms of the plaquette angle in the fol- 
lowing way: 

= -/3^cos6lp -7^cos26lp . (4.12) 
p p 

We use for the simulations the conventional AD hypercube with toroidal 
boundary conditions (hypertorus), and, for comparison we also consider 
the surface of a 5D cube which is topologically equivalent to a AD sphere. 
Contrarily to the hypertorus, the surface of a 51? cube is not homogeneous. 
There are a number of sites which do not have the maximum connectivity, 
namely 8 neighbors. Due to this fact, uncontrolled finite size effects are 
expected to turn up. Their influence can be somehow alleviated by the 
introduction of appropriate weight factors in those inhomogeneous sites 
[ p^ . Since the topology remains unchanged and we do not expect the 
order of the phase transition to be affected by the rounding, we do not use 
weight factors. 

As a notational remark, we label with L and N the side length for the 
7YT topology and for the HS topology respectively. 

Next we define the energies associated to each term in the action 

Ep^^CE^^ose^) , (4.13) 

^^p p 

E2p^^{J2cos2ep) , (4.14) 
p 

where Np stands for the number of plaquettes. On the hyper-torus 

this number is simply proportional to the volume, namely the forward 
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plaquettes are A^p = QL^. On the sphere the number of plaquettes has 
a less simple expression, and can be computed as a function of A'' as 
Np = 60{N - ly + 20{N - 1)2. In this case, the system is not homo- 
geneous and iVp is not proportional to the number of points on the four 
dimensional surface, which is — (iV — 2)^, some points having a number 
of surrounding plaquettes less than the possible maximum 12, as opposed 
to what happens on the torus . In order to allow comparison, we define 
LeS = {-^y^'^t in such a way that a hypertorus of i = Leff has the same 
number of plaquettes as the corresponding hypersphere. 

The specific heat and the Binder cumulant are useful quantities to mon- 
itorize the properties of a phase transition. At the critical point, they are 
known to posses different thermodynamical limits depending on whether 
the transition is first order or higher order, and hence, their behavior with 
increasing lattice size can give some clues as to the order of the phase tran- 
sition. We have studied the Finite Size Scaling (FSS) of these two energy 
cumulants. 

The specific heat is defined for both energies as: 

Cv = = N^iiE^) - {Epf ) , (4.15) 

C2V = ^E^^ = Np{{El) - {E^pf) . (4.16) 

As we have observed a high correlation between both energies, the re- 
sults being qualitatively the same, we shall only report on the observable 
defined for the plaquette energy E^. 

In a second order phase transition, scaling theory predicts the specific 
heat maximum to diverge as U^l" . If the transition is first order it is 
expected to diverge like the volume L'^ (or strictly like the number of pla- 
quettes iVp) reflecting that the maximum of the energy fluctuation has the 
size of the volume. This is expected to hold in the asymptotic region of the 
transition L ^c- In the transient region, namely L < ^c, the specific heat 
is expected to grow more slowly than the volume. 

We shall use the number of plaquettes iVp to study the specific heat 
behavior in order to have a single parameter for both the and the W<S 
topologies. This is of course equivalent to doing the discussion as a function 
of L on the UT and Leff on the US. We have in this case ~ (ATp)"/'*". 

We have also studied the behavior of the fourth cumulant of the energy 



{ED 
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When the energy distribution describing the system is gaussian Vl — > 
2/3 in the thermodynamic Hmit. This is the case for a second order phase 
transition. If the transition is first order, far from the critical couphng /3c, 
Vl also tends to 2/3, reflecting the gaussianity of the energy distribution. 
However, at /3c the distribution can be described by two gaussians cen- 
tered about the energy of each metastable state Ei and E2 , and hence this 
quantity has the non-trivial thermodynamic limit: 

Another interesting quantity measurable from the density of states is the 
distribution of the partition function zeroes, or Lee- Yang zeroes [|^. To 
clarify the role of the partition function zeroes let us first stress a well known 
fact: there are no phase transitions on a finite volume, phase transitions 
arise in the thermodynamical limit and are signaled by non-analyticity in 
the free energy: 

F(/3,y) = -^logZ(/3) . (4.19) 

Z{(3) is a linear combination of exponentials, and hence an analytical 
function. This implies that the free energy can be singular only where 
Z{(3) = 0. However if the coupling /3 is real, and the volume is finite, that 
linear combination is a sum of positive terms and hence, the zeroes of Z{P) 
are located in the complex plane of the coupling /3. The onset of the phase 
transition in the limit of infinite volume is signaled by a clustering of the 
zeroes on the real axis at /3c. 

Here follows the description of our procedure, for a general description 
of the method see . 

As mentioned previously we work only with Ep, and then we consider 
only the spectral density method p6| for this variable. Since we work at 
fixed 7 we consider only the /3 coupling derivatives. 

From the MC simulation at /3 we obtain an approximation to the density 
of states which allows us to compute the normalized energy distribution 
Pp{E). The energy distribution can be expressed as 

Pfi{E) = ^W{E)e-P'' . (4.20) 

We use the standard reweighting technique to obtain from the distribu- 
tion measured at /3 the distribution at another coupling oj, which is complex 
in the more general case w = 77 -|- z^. The standard reweighting formula is: 
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using ( 4.2Cl| ) and the normalization condition PpiE) = 1 one obtains 



^J2PpiE)e-^'^-'^^^ . (4.22) 



We can factorize the contributions from the real and imaginary part: 



Z{(3) 



J2 Pp{E){cosiE^) + i sm{E^))e-<-'^-^^^ . (4.23) 



This is the standard reweighting formula extended to the complex pa- 
rameter space of couplings. As a first observation we have a pure oscillating 
factor due to Im{uj) ^ 0. Since E is 0{V) this is a rapidly oscillating func- 
tion which makes it impossible to locate zeroes with large imaginary part. 

The real part of the coupling contributes to the well known exponential 
damping, e~^'^~^'^^ , which is telling us that we can trust the extrapolation 
only in a small neighborhood /3 ± ry. 

Since Z(Rbuj) has no zeroes in a finite volume, an easy way to locate 
the zeroes numerically is looking at the minima of the function |G(Li;)p, 
where: 



^(")-z(R^- t:m^)^^^^ ■ ^ ^ ^ 

The function to minimize is: 



_ {Y^^Pp{E)e-^^-^)'^cos{Ei)f + {Y^^Pp{E)e-^^-0)^sHEOf 

(4.25) 

One interesting property of the partition function zeroes concerns the 
estimation of critical exponents. Denoting by loq the coupling where the 
first zero is located, the distance to the real axis scales with the v exponent: 

Im{uJo) - L-^'" . (4.26) 

We shall use this property to compute the effective v exponent. 

From Pe{P)l, we can measure the free energy gap, AF(L), which is 
the difference between the minima and the local maximum of the free en- 
ergy |22]. We use the spectral density method to get, from the measured 
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Figure 4.1: Phase diagram of U(l) in the extended (/3,7) parameter space. The 
dotted lines have been taken from WM. The symbols correspond to our simulations 
on the torus. The crosses correspond to the peak of Cv{L = 16), the errors are 
not visible in this scale. The diamonds have been obtained with hysteresis cycles 
in L = 8. 

histograms, a new histogram where both peaks have equal height. We take 
the logarithm of those histograms and measure the energy gap. A growing 
AF{L) in the asymptotic region of the transition implies a first order phase 
transition. An increase proportional to L''"^ is expected For the tran- 
sition to be second order, AF{L) must stay constant with increasing lattice 
sizes. 

4.2.1 Schwinger-Dyson Equations 

As a further check we have implemented the Schwinger-Dyson equations 
(SDE) 1^ which allow one to recover the simulated couplings from the 
Montecarlo data. 

Let A{0) be an operator with null expectation value: 

{A{0)) = Z-^ J [rf6i]A(6i)e-^[^l EE . (4.27) 

Derivating with respect to 9 this equation trivially yields 
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which is the equation of movement of the operator A(9) or Schwinger- 
Dyson equation. When the action depends on several couphngs the equa- 



tion (4.2g) can be expressed as: 



i 

This equation relates the values of the couplings with the expectation 
values we measure from the MC simulation. We need as many independent 
equations as couplings we have to determine in the action. In our case, in 
order to measure both /3 and 7, we need two operators with null expectation 
value in order to have two independent tests. At each lattice site n and for 
every direction /i we consider the operators: 

A{e) = sinOp = sin(0„,^ - Ostap) , (4.30) 
B{e) = sin 2ep - sin 2(^^„,^ - estap) , (4.31) 

where 9stap is the st aple of the link labeled by {n,fi). 
Applying equation ( 4.29[ ) to those operators we get: 



cos(6l„,^ - Ostap)) = /3(^ sin(6l„,,i - 9stap) ^ sin{9„,p - Ostap)) + 
p p p 

+27{^ sin(6'„,^ - Ostap) sin 2(6»„,^ - Ostap)) (4.32) 
p p 

2cOs2(e„,^ - Ostap)) = Sin2(6'„,^ - Ostap) ^ Sin(6l„,,, - Ostap)) + 

p p p 

-|-27{y^ sin 2{0„,^, - Ostap) ^sin2(e„.^-e,t<,j,)|4.33) 
p p 

On the hypertorus J2p means the sum over the plaquettes in positive and 
negative directions (±/x) bordering the link (12 plaquettes). On the sphere 
one has to be careful, since not all the links have 12 surrounding plaquettes, 
and the sum has to be understood as extended to the existing plaquettes. 

These equations hold for all n in such a way that we can sum up the 
equations for e very single link and average over the number of links. We 



quote in Table 4.1 for the tests done at 7 = —0.2 using both topologies. 
They show a perfect agreement between the simulated couplings and the 
ones recovered from the MC expectation values. 
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HT topology 


TiS topology 


L 


/5s im 


(,fl,7)sD 


N 


/3sim 


(/3,7)sD 


6 


1.1460 


1.1466(44),-0.2009(26) 


6 


1.1587 


1.1592(14),-0.1998(12) 


8 


1.1535 


1.1527(22) -0.1997(15) 


8 


1.1600 


1.1599(12) -0.1999(8) 


12 


1.1582 


1.1581(12),-(). 1999(8) 


10 


1.1602 


l.lG02(7),-0. 2001(9) 



Table 4.1: Couplings obtained from the MC simulations at 7=-0.2 using the 
Schwinger-Dyson equations. 

4.3 Numerical Simulation 



7 = +0.2 


N 




/5sim 


r 


iVr 


/?*(i) 


Im(a;()) 


4 


5.383 


0.8910 


70 


29142 


0.8906(2) 


0.0169(1) 


5 


7.150 


0.8855 


210 


8000 


0.8855(2) 


0.00604(6) 


6 


8.921 


0.8834 


520 


3800 


0.88330(5) 


0.00244(6) 


7 


10.694 


0.8818 


790 


2300 


0.88182(3) 


0.00090(3) 


8 


12.469 


0.88095 


930 


1930 


0.88097(3) 


0.00040(2) 


7 = 


N 




/3sim 


r 


iVr 


/3*(i) 


Im(ci;o) 


6 


8.921 


1.0128 


240 


8300 


1.01340(5) 


0.00616(8) 


8 


12.469 


1.0125 


400 


2300 


1.0127(2) 


0.00196(4) 


10 


16.021 


1.0120 


780 


1200 


1.01212(3) 


0.00075(3) 


12 


19.574 


1.0119 


850 


1150 


1.01194(2) 


0.00031(2) 


14 


23.123 


1.0117 


920 


1900 


1.01168(2) 


0.00013(1) 


7 = -0.2 


N 




/5sim 


r 


iVr 


/3*(L) 


Im(a;o) 


6 


8.921 


1.1587 


160 


2000 


1.1587(4) 


0.0107(3) 


7 


10.694 


1.1600 


320 


2800 


1.1597(1) 


0.00550(7) 


8 


12.469 


1.1597 


510 


1000 


1.1603(2) 


0.00321(4) 


10 


16.021 


1.1602 


680 


1500 


1.1604(2) 


0.00171(2) 


12 


19.574 


1.1604 


820 


1200 


1.1602(1) 


0.00083(1) 


14 


23.123 


1.1605 


900 


1100 


1.16048(5) 


0.00044(1) 


16 


26.684 


1.1604 


1150 


1200 


1.16038(2) 


0.00023(1) 



Table 4.2: Statistics of the data obtained for the HS topology 
Most of the work has been done by simulating the subgroup Z(1024) c 
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7 = -0.1 


L 




T 






7m (wo) 


6 


1.0720 


350 


1900 


1.0716(2) 


0.0097(1) 


8 


1.0784 


640 


1400 


1.0786(2) 


1.1539(2) 


12 


1.0820 


820 


750 


1.0818(1) 


0.00114(2) 


16 


1.08278 


930 


900 


1.0827(1) 


0.00040(2) 


20 


1.0833 


1150 


1100 


1.0833(1) 


0.00020(1) 


7 = -0.2 


L 


/5s im 


r 


iVr 




Im[u)a) 


6 


1.1460 


380 


2000 


1.1452(2) 


0.0115(1) 


8 


1.1535 


620 


1900 


1.1539(2) 


0.00506(6) 


12 


1.1582 


840 


1200 


1.1582(2) 


0.00152(3) 


16 


1.15935 


920 


900 


1.1593(1) 


0.00060(2) 


20 


1.1599 


1150 


900 


1.1599(1) 


0.00028(1) 


7 = -0.3 


L 


/5sim 


T 






Jm(a)o) 


6 


1.2255 


340 


2000 


1.2237(4) 


0.0138(1) 


8 


1.2344 


560 


1800 


1.2340(1) 


0.00623(4) 


12 


1.2395 


770 


1000 


1.2395(2) 


0.00190(3) 


16 


1.2410 


900 


900 


1.2410(1) 


0.00084(2) 


20 


1.2416 


1100 


1200 


1.24156(5) 


0.00041(2) 


24 


1.2417 


1200 


1100 


1.24162(5) 


0.00022(1) 


7 = -0.4 


L 




r 




P*{L) 


Jm(o;o) 


6 


1.3090 


400 


1800 


1.3082(4) 


0.0152(1) 


8 


1.3192 


600 


1600 


1.3194(3) 


0.00704(5) 


12 


1.3258 


710 


1400 


1.3259(1) 


0.00238(3) 


16 


1.32775 


840 


900 


1.3278(1) 


0.00108(2) 


20 


1.3285 


930 


1600 


1.3284(1) 


0.00054(2) 


24 


1.3286 


1150 


1500 


1.3286(1) 


0.00029(1) 



Table 4.3; Statistics of the data obtained in the T-CT topology. 



f7(l), since the phase transition associated to the discrete group hes safely 
far away. The overrelax algorithm can be applied only in the simulations 
with 7 = so that the gain in statistical quality due to the overrelax effect 
is limited to the Wilson action, whore wo have simulated both, the full 
group U{1) and the discrete one for the sake of comparison. For 7 ^ we 
have simulated the discrete group since so the simulation is considerably 
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speeded up. 

For every lattice size we consider, we perform trial runs to locate the 
coupling P* (L) where the specific heat shows a peak. We use the standard 
reweighting techniques to extrapolate in a neighborhood of the simulated 
coupling. Once the peak is located within an error in the fourth digit of the 
coupling, we perform an intensive simulation there to get Pe {(3* ) l ■ The 
statistics performed at these pseudo-critical couplings are reported in Table 



4.2 and Table 4.3 for the TiS and the Ti.T topology respectively. 



Typically we measure the energies every 10 MC sweeps in order to 



construct the energy histogram (4.20) From this energy distribution we 



obtain the cumulants of the energy we are interested in, and the critical 
exponents using Finite Size Scaling techniques. We remark that one of the 
main sources of systematic error when measuring critical exponents, is the 
indetermination in the coupling [3* (L) where to measure. 

The simulations have been done in the RTNN machine consisting of 32 
Pentium Pro 200MHz processors. The total CPU time used is the equivalent 
of 6 Pentium Pro years. 

We have updated using a standard Metropolis algorithm with 2 hits. 
The acceptance has always been between 65% and 75%. 

In order to consider the statistical quality of the simulation, following 
[E7| we define the unnormalized autocorrelation function for the energy 



N-t 



1=1 

as well as the normalized one 

Mt) ^ §1 . (4.35) 

The integrated autocorrelation time for the energy, r™*, can be mea- 
sured using the window method 

T™*(t)-^ + E'^(*') ' (4-36) 

t'=i 

for large enough t, which is in practice selected self-consistently. We use 
t in the range 5r*"*, lOr™*, and we check that the obtained r™* remains 
stable as the window in t is increased. 

We have always started from hot and cold configurations in order to 
make sure that the system does not remain in a long living metastable 
state, which could be interpreted as a Dirac sheet. The results coming 
from both types of starts have always been indistinguishable. 
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4.4 The phase transition hne /3c (7) 



We have studied the deconfinement-confinement phase transition hne at 



several values of 7 (Figure 4.1) 



In the region of negative 7, apart from the deconfinement transition 
hne, there is another transition hne provoked by the competing interaction 
between the couphngs {(3,j). In the hmit 7 = —00 the model is not dy- 
namical since for all finite /3, cos 6'p = 0. So we expect this line to end at 
the corner (7 = — 00,/?= +00). We have not performed a deep study of 
this transition line, however the simulations at L = 8 revealed double-peak 
structures pointing to a first order character. The limit 7 = -|-oo is equiv- 
alent to a Z{2) theory and the critical point can be calculated exactly by 
self-duality. 

We have studied carefully the region between lines A and B finding no 
signatures of the existence of additional lines. 

We have focused on the transition line A, at several values of 7, on the 
hyper-torus and on the hyper-sphere. 

The structure of this section is the following: 

First, we describe the results at 7 = +0.2 on the spherical topology. 
Our purpose is to observe how does the TiS topology behave in a non con- 
troversial region, and comparing with the known results on the hypertorus. 

Second, we study the point 7 = (Wilson action) on the spherical 
lattice. This coupling has been recently claimed to be the starting point of 
a critical line (infinite correlation length at the critical point) which extends 
in the range 7 < 0, with an associated critical exponent v ^ 0.37 [Q. Since 
this assertion relies on the absence of two-state signals on spherical lattices 
up to iV = 10, we shall check whether or not double peak structures set in 
when larger spheres are considered. 

Third, we go to the region of negative 7 values. We are aware of no 
previous systematic study on the hypertorus in this region, so we run sim- 
ulations up to L = 24 at 7 = -0.1,-0.2,-0.3,-0.4. Motivated by the 
results on the hyper-torus we have just run a single 7 negative value on the 
hyper-sphere. For the sake of comparison with we choose this value to 
be 7 = -0.2. 

Finally, we discuss the Finite Size Scaling behavior exhibited by both 
topologies. 



4.4.1 Results at 7 = -1-0.2 on the spherical topology 

As we have pointed out, there is a general agreement on considering the 
phase transition first order in the region of positive 7. 
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Figure 4.2: MC evolution of Bp at 7 = +0.2 on the HS topology for iV = 6 
(lower window), A?^ = 7 (middle) and N = 8 (top). 
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Figure 4.3: Ep distribution at 7 = +0.2 on the HS topology. The small window 
is the distribution we obtained on the TiT topology in L = 4 at /3 = 0.8595 
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Figure 4.4: Binder cumulant at 7 = +0.2 on the TiS topology in N = 4, 5, 6, 7, 8. 
The cross in the upper corner signals the second order value 2/3. 
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Figure 4.5: Specific heat maximum, and Ferrenberg-Swendsen extrapolation 

(solid line) at 7 = +0.2 on the 7^5 topology. The small window represents 
C™^^(Afp) The dotted hne corresponds to the slope expected in a first order 
phase transition. 
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Figure 4.6: Ep distribution at 7 = +0.1 on the TtS topology. 



The lattice sizes used range from = 4 to iV = 8, which correspond to 
Lc« ~ 5 and LcS ^ 12 respectively see (Tabic 4.2). 

In Figure 4.2, the MC evolution for = 6,7,8 is shown. We remark 
that no multicanonical update is needed to obtain a very high rate of flip- 
flops up to iV = 8. However, on the torus, for L > 6 the probabihty 
of tunneling between both metastable states is so tiny that a reasonable 
rate of flip-flops is not accessible to ordinary algorithms. Probably the 
inhomogeneity of the sphere is increasing the number of configurations with 
energies which correspond neither to the confined nor to the deconfined 
phase, but in between. These configurations make the free energy gap to 
decrease and hence the tunnehng is easier on the HS topology. In short, 
the sites without maximum connectivity act as catalysts of the tunneling. 

In Figure 4.3 the energy distributions are plotted. The distribution 
at TV = 6 is distinctly non-gaussian, but a blatant two-peak structure is 
observed only from N ~ 7 {Left ^ 11) on. So, when comparing with the 
result at 1/ = 4 on the torus, (see small window in Figure |4.3| ) a first 
observation is that the onset of a two-state signal is particularly spoiled by 
the HS topology, at least in AD pure compact U(l) gauge theory. 

Concerning the latent heat, we remark its stability aheady at iV = 7, 
or, if anything, its increase from N = 7 to N = 8. One would even say 
that the same happens between the positions of the would-be two states in 
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iV = 6 and the position observed in iV = 7. The behavior of the Binder 
cumulant reflects this fact (see Figure iA). We observe a rapid growth of 
y^nin ^^j. gj-^g^^ lattices sizes, apparently towards 2/3 (second order value). 
This growth stops when the splitting of the two peaks is observed. The 
splitting of the two peaks is reflected by smaller values in V^"". We remark 
that the errors quoted for V2"'" are calculated taking into account the in- 
determination in the value of the minimum, but not the displacement in 
the position of the coupling where the minimum appears. 

In our opinion, mainly two reasons can give account of this behavior. 
The first one comes from general grounds: at the very asymptotic region 
of a first order phase transition, the energy jump gets sharper and sharper, 
and a slightly increase of the latent heat could be expected. The second 
one would be the increasing restoration of homogeneity in the hypersphere 
when increasing the lattice size. The latent heat observed in small N might 
be affected by the inhomogeneity of the hyper-sphere. 

At this point we cannot give a single reason for this to happen. We 
postpone a stronger conclusion to the section devoted to the Finite Size 
Scaling discussion. 

The value of the latent heat when obtained with a cubic spline fit to the 
peaks in = 8 is Ciat = 0.064(2). Results obtained with mixed hot-cold 
starts in A = 9 seems to give an energy jump around 0.067. However the 
reliability of this method is very limited and we do not dare to extract 
strong conclusions from it. Taking into account that he cubic spline at 
N — 7 gives 0.053(2) a possible scenario would be a slowly increasing latent 
heat towards its asymptotic value Ciat(oo). 

The peak of the speciflc heat for the different lattice sizes is displayed in 
Figure Lb. The continuous line represents the FS extrapolation. The plot 
of the peak value, C^^^{Np) as a function of the plaquette number reveals 
a linear relationship, and hence a flrst order character. 

We have also run simulations at 7 = -fO.l in order to get an estima- 
tion for the latent heat. The results, being qualitatively similar to those 
encountered at 7 = -t -0.2 are not reported in detail. We plot the energy 
distribution in Figure i.t. 



4.4.2 Results at 7 = on the spherical topology 

Let us flrst situate the status of the studies using the toroidal topology 
with the Wilson action. As we have pointed out in the introduction, the 
transition was believed to be continuous till the simulation of a L = 6 
lattice revealed the flrst signs of the existence of two metastable states 
Numerical simulations up to L = 10 showed a rapid decrease of the latent 
heat when increasing the lattice size, suggesting that the two metastable 
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states might superimpose in the hmit L — s- oo However, in Azcoiti 
et al. suggest that the latent heat starts to stabihze at L ^ 12, 14. 

When one tries to simulate large lattices the tunneling becomes scarce 
and the use of multicanonical simulations for those lattice sizes seems to 
be in order. However, in spite of its usefulness in spin models, as far as we 
know, there are no multicanonical simulation showing a substancial flip-flop 
rate improvement for this model. 

Lattice sizes up to L = 16 have also been studied using multihistogram- 
ming techniques |^ and RG approaches ||l^, 11 . The results support the 
idea of a quasi stable latent heat for L > 12. 

Topological considerations stressed already in the introduction led some 
authors to use lattices homotopic to the sphere. In the case of the Wilson 
action (7 = 0) the two-state signal is absent up to iV = 10 |l^. Following 
this, we have simulated on the spherical topology in lattices ranging from 

= 6 to iV = 14, finding that the two-state signal sets in fro m N — 12 on. 

The MC evolution for N = 10, 12, 14 is plotted in Figure |j[ For the 
TV = 14 lattice we have four independent runs, signaled by the dashed lines 
in the figure, all of them giving the same predictions. 

In Figure 4.8 the distribution of Ep is plotted in lattices ranging from 
iV = 6 to iV = 14. We observe that a two-peak structure is revealed first 
time by the histogram in iV = 12, which has a L^s ~ 19. On the toroidal 
topology the equ ivale nt signal is observed already at L = 6 (see small 
window in Figure |4.§| ) . 

Having in mind the results of the previous section at 7 = +0.2, it seems 
that the minimum L^s required to observe a two-peak structure in the 
spherical topology is around three times the minimum lattice size L needed 
in the torus to observe two peaks. 

The lattice sizes used in [|8j at 7 = ranged from = 4 to = 
10, so it is not surprising that they did not see any two-peak structure. 
Our results are in that sense compatible with theirs, though ours show a 
faster divergence for C^'^^{N) with increasing A^ than the one observed in 
[Q in the lattice sizes we share (N=6,8,10). However we are particularly 
confident on this point because our simulation has been performed closer 
to the peak of the specific heat P*{L), and so we expect the Ferrenberg- 
Swendsen extrapolation to be more precise. 

For the sake of comparison with the full group U(l), simulations with 
the Wilson action have been performed at 7 = 0. At this coupling an 
implementation of overrelax is possible, and the global decorrelating effect 
should manifest itself in a better statistical quality. 

In Figure 4.9 we show the histograms of Ep from A^ = 6toA^ = 12 sim- 
ulating the full group. The results are fully compatible with those obtained 
from simulations with Z(1024). 
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Figure 4.7: MC evolution of Sp at 7 = on the HS topology for TV = 10 (lower 
window), A*' = 12 (middle) and A'' = 14 (top). The vertical dashed lines separate 
different runs 



This test being performed, we go back to the description of the results 
obtained for the discrete group. 



In Figure 4.10 we plot Cv(TV) for different A^'s, together with the Ferrenberg 
Swendsen extrapolation. The small window is a log — log plot of C^'^^{N) 
as a function of A^p. A linear behavior is observed from A^ = 10 on, which 
is even faster than linear when A^ = 14 is taken into account. This fact 
by itself implies the first order character of the transition since it means 
that the maximum in the energy fluctuation has the size of the volume. As 
a further check we have measured the v exponent from the scaling of the 
Fisher zeroes (see Table 4.2). 

In Figure 4.11 we plot |G(ci;)p for several lattice sizes on the spheri- 
cal topology at 7 = 0. The different curves stand for the different oj we 
extrapolate. For small lattice sizes Im{ujQ) is larger, and the damping is 
more severe than for the larger lattices since in the later the imaginary part 
contributes with a faster oscillating function. Actually, we observe that for 
A^ = 14 even a second minimum could be measured before the signal is 
damped, while for N = 10 one can measure accurately only the first one. 

The results for the effective ly are plotted in Figure 4.12 lower win- 



dow. As could be expected from the behavior of the specific heat, for small 
lattices the effective exponent is somewhat larger than 0.25. It gets com- 
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Figure 4.8: Ep distribution at 7 = on the HS topology. The small window 
corresponds to our simulation in the HT topology in L=6 at /3 = 1.0020. 




Figure 4.9: Ep distribution at 7 = on the HS topology simulating the full 
group U(l). 
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Figure 4.10: Specific heat maximum, and Ferrenberg-Swendsen extrapolation 
(solid line) at 7 = on the TiS topology. The small window represents C™^''(A'^p) 
The dotted line corresponds to the slope expected in a first order phase transition. 



patible with the first order value from = 10 on. In the largest lattices 
the distance between the two peaks slightly increases, and we measure a 
VcS slightly smaller than the first order value. As a < 0.25 is strictly 
impossible we expect this to be a transient effect due to finite size effects 
associated to the observed splitting up of the two peaks. 

Concerning the latent heat, we observe a behavior completely analogous 
to the one observed at 7 = +0.2. We do observe a two peak structure quite 
stable when comparing iV = 12 and N = 14, but it slightly increases in 
N = 1 4. T he plot of the Binder cumulant reflects again this fact (see 
Figure 4.15 ). The fast growth of V£^™ towards 2/3 is preempted by the 
onset of double peaked histograms from = 12 on, and it even decreases 
in TV = 14. 

The cubic spline fit in = 14 at the peaks gives a latent heat Cut ~ 
0.018(2). which is compatible with the results suggested by extrapolating 
to infinite volume the values obtained on the torus up to L = 14. Our 
results are hence supporting the conjecture stressed in about a quasi 
stabilization of the latent heat on the torus from L = 12 on. 
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Figure 4.11: Plot of ( [4.25| ) for iV = 10 (lower window), TV = 12 (middle) and 
N =14 (top) at 7 = in the HS topology. 



7 = -0.2 



4.0-10" 
1/Np 



7=0 



4.0-10" 
1/Np 



Figure 4.12: Effective exponent v a,t ^ = (lower window), and at 7 = —0.2 
(upper window) on the spherical topology. In the smaller lattices an i^cff ~ 1/3 
is observed which becomes 1/d when large enough lattices are considered. 
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Figure 4.13; Binder cumulant at 7 = (Wilson action) on the HS topology in 
N = 7,8, 10, 12, 14. The cross in the upper corner signals the second order value 
2/3. 



Figure 4.14: MC evolutions at 7 = 



-0.2 in L=12,16,20 on the HT topology. 
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Figure 4.15: Ep distributions at the diflcrcnt negative 7 on tiie TCT topology. 
We have run up to L = 20 at 7 = -0.1, -0.2 and up to L = 24 at 7 = -0.3, -0.4. 
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Figure 4.16: Free energy gap AF{L) (lower window) and latent heat (upper 
window) for the different negative 7 in the TiT topology. 
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Figure 4.17: Effective v exponent for the different negative 7 in the HT topology. 
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Figure 4.18: C^'^'^{L) as a function of A'p for the different negative 7 values on 
the TiT topology. We have used log — log scale for the clarity of the graphic sake. 
The dotted line corresponds to the first order slope. 
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4.4.3 Results for 7 < 



Toroidal topology 

We have performed a systematic study of the transition at several negative 
values of 7 on the toroidal topology. 

We find that the two-state signal persists for all 7 values we consider. 

The MC evolution at 7 = —0.2 is plotted for several lattices sizes in 



figure 4.14 . As a general observation the more negative is 7, the fiip rate 
becomes scarce for increasing lattice size. As an example, at 7 = —0.4 the 
flip rate in L = 20 is comparable to the one observed in L = 16 at 7 = —0.2. 

The energy histograms reveal an increasing weakness of the transition 
when going to more negative 7 values (see Figure 4.15| ). A double peak 



structure is clearly visible at7 = — 0.1inL = 8, while at 7 = —0.4 one has 
to go to L = 12 to observe an equivalent signal. 

In what concerns the behavior of the free energy gap AF(L), it grows 



for all investigated lattice sizes at all negative 7 values (see Figure 4.16 
lower window) signaling the first order character of the transition. The 
value of L at which AF{L) starts growing is certainly larger as the value of 
7 is more negative. This is another test of the increasing weakness of the 
transition as 7 gets more negative. 

The statistics performed on the torus are reported in Table |4.3| . We also 
quote for the diS'erent negative 7, the value of P*{L) and the position of 
partition function zero closest to the real axis. We have computed from the 
imaginary part of the zeroes the effective v exponent between consecutive 



lattice sizes following (4.26) 



In Figure 4.17 we plot for the different 7 values the we measure. In 
all cases a v^a ~ 1/3 is observed for small lattice sizes, which gets closer 
to 0.25 when the lattice is large enough. From this figure the trend of v^s 
seems rather clear towards the first order value. 

From the energy distributions we measure the latent heat through a 
cubic spline at the peaks. Taking into account the value v^s = 0.25 we 
measure, we plotted it as a function of the inverse of the volume 
which is also the expected behavior of the latent heat when the transition 
is first order. The latent heat can be extrapolated to a value which is safely 



far from zero (see Figure |4T6| upper window) 



In Figure 4.18 we plot C^'^^{L) as a function of A'p. As could be ex- 
pected from the behavior of the effective exponent the maximum of 
the specific heat for small lattices diverges slower than the volume. For 
the smaller lattices the effective exponent is a/iy ^ 1.4. This value in- 
creases monotonically with the lattice size. In the largest ones we observe 
a/v ~ 3.5. 
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Figure 4.19: MC evolution of i5p in iV = 12 (lower window), iV = 14 (middle 
window) and A'^ = 16 (upper window) at 7 = —0.2 on the Ti.S topology. The two 
different runs in A*' = 16 are separated by a vertical dashed line. 



Spherical topology 

On the toroidal topology we have found that the minimum lattice size 
required to observe a two state signal is obtained through an apropriate 
combination (7, Lmin), with increasing Lmin for increasingly negative 7, the 
behavior being qualitatively similar for all the 7 values we have investigated. 
In view of this, we have studied a single 7 < value on the spherical lattice 
to check if the two state signal is absent on this topology. For the sake of 
comparison with |l^, |l^ we choose this value to be 7 = —0.2. 

We have run simulations on spheres ranging from N = 6 {Left ~ 8) to 
iV = 16 (Loff 26). The MC evolution is plotted in Figure |l9. 

The distribution in = 14 is distinctly non gaussian and the splitting 
of the peaks occurs in = 16 (see Figure 4.20). We remark that the 
simulations in N — 16 are extremely expensive in CPU. In order to alleviate 
thermalization we have parallelized the code using shared memory in two 
PPro processors. In this lattice we have run two independent simulations 
starting from different configurations the results being fully compatible. 
The behavior of the maximum of the specific heat is shown in Figure 
We observe the same trend than in the previous values of 7 in the 
sphere, i.e., an increasingly fast divergence of the specific heat with the 
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Figure 4.20: Ep distributions at 7 = —0.2 on the HS topology. 




Figure 4.21; Specific heat maximum, and Ferrenberg-Swendsen extrapolation 
(solid line) at 7 = —0.2. The small window represents C™^''(iVp) The dotted line 
correspond to the slope expected in a first order phase transition. 
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Figure 4.22: Binder cumulant at 7 = —0.2 on the TiS topology in 
A'' = 10, 12, 14, 16. The cross in the upper corner signals the second order value 
2/3. 



lattice size, showing an effective a/i^ ~ 4 already when the histogram width 
becomes constant with increasing lattice size. 

At this value of 7 we have a worse estimation for the latent heat to be, 
as the two peaks have not split enough to allow an accurate measurement. 

Agai n the behavior of the Binder cumulant is very significant (see Fig- 
ure |4.22| . The value y^°'" shows a very fast trend towards 2 /3 in the small 
lattices. When increasing the lattice size the rate gets slower, and finally 
the value in = 16 is compatible with the value in A^ = 14 preempting the 
extrapolation to 2/3. Unfortunately lattices larger than A^ = 16 are unac- 
cessible to our computers nowadays, and we cannot observe a decreasing 
V^™ as we did for the other 7 values. However from the behavior between 
A^ = 14 and A^ = 16 an increasing V^™ for larger spheres seems to be very 
unlikely. 

Concerning the effective critical exp onen t v we have measured the po- 
sition of the first Fisher zero (see Table 12 ) and compute i/eS (see Figure 



4.12). Again the first order value l/d is reached in the largest lattices. 
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4.4.4 Toroidal versus Spherical topology 



Due to the fact that on the sphere there are a number of points with less 
than the maximum connectivity, violations to standard FSS in the form of 
uncontrolled finite size effects are expected to appear on this topology. 

We have observed all along this work that working on the spherical 
topology retards the onset of two-state signals. In terms of Leff, some 
qualitative prediction arising from our results would be that at fixed 7, the 
minimum lattice size required to observe a double peak structure, Lmin, is 
around three times larger on the sphere than on the torus. This one being 
the showier difference, is not certainly the only one. 

When comparing C™'^^{Lcs) on both topologies one finds alway s smaller 
values on the spherical topology (see small window in Figure 4.23 ). We did 
run on L = 12 on the torus and TV = 8 on the sphere, which has an 
LcS ^ 12, sweeping an interval of /3 values including the phase transition. 
In Figure 4.23 we have plotted Ep in both topologies. In the region of low 
(3 the system is disordered, the entropy is higher and the system is not so 
sensitive to inhomogeneities of the lattice. However, in the region of high /3 
the energy is in general smaller than for the homogeneous lattice since the 
system tends to be ordered, and the influence of the sites with less than 
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Figure 4.23: Ep at j = —0.2 in Leg ~ 12. The dotted line corresponds to the 
toroidal topology, the dashed one to the sphere N — 8. The small window shows 
the difference in the specific heat between both cases. 
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maximum connectivity is more evident. This difference gets smaller when 
the lattice size is increased on the sphere, and hence the energy jump is 
larger. This feature explains the observed splitting of the two states on the 
spherical lattices when increasing lattice size. 

Altogether, it is not recommended to work on spherical lattices to check 
the existence or not of two-state signals. However, there are a number of 
facts that make this topology not so dissapointing in spite of those uncon- 
trolled finite size effects. The first one concerns the behavior of /3* (L) and 
the second one the measure of the latent heat. 



In Figure 4.24 we compare (3*{L) on the torus and on the sphere. A first 
observation is that despite they not having the same values in finite lattices, 
both curves get closer when increasing the lattice size. This supports the 
idea of a common thermodynamic limit for both topologies. On the other 
hand, the shift in the apparent critical coupling with L is much less dramatic 
on the sphere than on the torus. It seems that finite size corrections to j3c 
are smaller in the case of the spherical topology. This behavior has also 
been observed in the Ising model on spherical lattices pst . 

The behavior of the latent heat is very significant. In Figure 4.25 we 
compare the latent heat in finite volumes for both topologies. On the 
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Figure 4.25: Latent heat for the different lattice sizes on the torus and on the 
sphere at 7 = (lower window) and at 7 = —0.2 (upper window). The dotted 
line corresponds to Ciat(oo) obtained by extrapolating the values on the torus 
through a linear fit as function of 1/V . The values for 7 = on the torus have 
been taken from 



Wilson line, at 7 = the asymptotic value of the latent heat is obtained in 
the spherical topology. The inhomogeneity of the sphere has, paradoxically, 
helped us to run lattices with L^s ~ 24 without having to worry about the 
tunneling rate. From this graph we quote Ciat(oo) — 0.018(1). In the upper 
window the same is plotted for 7 = —0.2. Unfortunately we have just a 
single lattice size on the sphere to measure the latent heat, however, from 
the behavior exhibited by this topology, an increase of the latent heat on 
spheres larger than = 16 seems rather likely. Altogether, at 7 = —0.2 a 
Ciat(oo) ^ 0.009 is plausible. 

In what concerns the latent heat, the spherical topology seems to afford 
an useful complement to the results obtained for the the toroidal lattice. In 
fact, on the Wilson line lattices larger than L = 16 cannot be studied due 
to the technical problem associated to the scarce tunneling. The spherical 
lattice alleviates this technical problem, and the results are supporting the 
value of the extrapolated latent heat as a function of the inverse of the 
volume, from the data on toroidal lattices up to L = 16. One could consider 
the possibility of simulating on spherical lattices at negative 7 to solve the 
tunneling problem, but in our opinion, the price to pay is too high. 
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4.5 Conclusions 



The first order character of the deconfinement transition in pure U(l) has 
been proved, up to the hmits of a rather reasonable numerical evidence, in 
the interval 7 G [+0.2, -0.4]. 

In 7 = +0.2 we have been able to stabilize the latent heat. We are aware 
of no simulation on the torus showing a stable latent heat due to the scarce 
tunneling in lattices larger than L = 6. The spherical topology has helped 
to solve this problem. However, probably any lattice with inhomogeneities 
would produce the same catalysing effects. 

In 7 = we have obeserved an increasing energy gap on the spherical 
lattice, and also, we have been able to measure the asymptotic value of 
the latent heat. We have proved the suggestions of several authors about a 
quasi stabilization on lattices larger than L ^ 12. The data on the spherical 
topology have been crucial to discard the possibility of a slowly vanishing 
latent heat. It follows that the discretization of pure compact U(l) LGT 
on the lattice using the Wilson action exhibits a first order phase transition 
with a latent heat in the thermodynamical limit Ciat ~ 0.018. 

On the toroidal topology things happen qualitatively in the same way 
than in 7 = up to 7 = —0.4. We have run on spherical lattices in 
7 — —0.2 looking for an argument similar to the one found in the Wilson 
case concerning the stabilization of the latent heat. The simulation in 
A'^ = 16 gives an estimation for the latent heat though rather imprecise 
because the splitting of the peaks is not good enough for the measurement 
to be accurate. However, in view of the behavior exhibited by the spherical 
topology concerning the trend of the two peaks on spherical lattices to 
separate, we are prone to consider the value measured in A'^ = 16 as a 
lower bound for the latent heat in 7 = —0.2. On the other hand, we 
would at present say that the latent heat extrapolated from the data on 
the torus up to L = 20 is rather accurate, in view of the behavior for 
the Wilson case. The possibility of running larger spheres surpasses our 
computer resources. A highly parallelized version of the code should be 
used to alleviate thermalization, which is a possibility we do not discard 
completely at medium term. 

As for the evolution of the latent heat along the transition line as a 
function of 7. the scenario that follows from our results is plotted in Figure 

Either proving or discarding the possibility of a TCP at some finite 
negative |7tcp| > 0.4 will be a very difficult task from the numerical point 
of view. An analytical argument would be welcome. A very tiny two-state 
signal, comparable with the one observed in L = 8 at 7 = —0.4, is observed 
on the torus in L = 16 at 7 = —0.8, which is the most negative value we 
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Figure 4.26: Extrapolated value for the latent heat as a function of 7. The 
values quoted for 7 < come from data on toroidal lattices, the ones for 7 > 
from data on spherical lattices. 

have run on the torus. Taking into account the factor 3 in L^s one would 
need a A'^ '--^ 30 sphere to observe a tiny double peak structure. Simulation 
with spherical lattices in this range of 7 have no sense, and nothing can 
be concluded from the absence of two state signals from such negative 7 
values. 

In view of these difficulties to stabilize the latent heat for very negative 
7 values, the only chance to discern the order of the phase transition is 
the study of effective critical exponents. Cumulants of the energy, such as 
the Binder cumulant, have been shown to behave like expected in the first 
order case only when the stable two-peak structure is almost setting in, and 
there, we do not need further evidences any longer. We could not expect 
additional information since by definition the Binder cumulant relies on the 
existence of stable latent heat in order to extrapolate to V[™='oo ^ 2/3. 

The advantage of studying the effective exponents (which is nothing 
but studying the evolution of the histograms width) is that we do not need 
a direct observation of latent heat to conclude that a transition is first 
order js^, . We have observed a VcS which evolves monotonically until it 
reaches the first order value, namely 0.25, in all cases. The statistics needed 
to observe a monotonous behavior, are order IOOOt at the coupling f3*{L), 
which has to be located with high precision (four digits in our experience) 
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in order to accurately measure effective exponents. However, one has to 
make every effort to observe tlie trend of the effective exponents, since it is 
the only chance to discern the order of such tricky transitions. 

In pure U(l), for the Wilson case, an exponent ~ 1/3 was widely 
predicted ^, ^, ^ at the beginning of the eighties when the lattice sizes 
where too small to reveal two peaks. That v was shown to become v ~ 0.29 
when simulating L = 14 0. We have measured v = 0.25 on spherical 
lattices and stated its first order character. Within the approximation of 
effective potentials it can be shown that along the transient region of a weak 
first order phase transition, everything goes like in a second order one with 
a thermal exponent a = 0.5 Together with Josephson law (a = 2 — i/c?) 
it implies v ~ 0.37. This statement has been checked in 2D Potts |^^, 3D 
and 4D 0{N) models |§, |o| 



|31|, ^ and in the 4D SU(2)-Higgs at T=0 
with fixed Higgs modulus [33|. Our results prove that pure compact U(l) 
theory behaves in the same way. 
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Chapter 5 

Modelos bidimensionales 
de flujo de particulas 
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5.1 Introduccion 



En los liltimos afios esta recibiendo una especial atencion la formulacion de 
modelos para describir por ejemplo el paso de mensajes a traves de redes de 
ordenadores , el trafico en grandes ciudades HI , la adsorcion de moleculas 
en un cristal conductores de iones rapidos ||^. 

Todos estos problemas tienen en comiin que pueden ser estudiados desde 
el punto de vista de la teon'a de Random Walks. Sin embargo, el ele- 
vado numero de grados de libertad hacc imposble extraer informacion rele- 
vante mediante calculos anali'ticos. La alternativa es el uso de simulaciones 
numericas. 

Desde este punto de vista uno podria pensar en formular modelos tan 
proximos a la realidad como sea posible, incluyendo todos los grados de 
libertad que se sea capaz de manejar. Sin embargo los modelos que resultan 
aplicando esta filosoffa son demasiado complicados de estudiar, no solo en 
lo que concierne a la descripcion de fenomenos relevantes, sino tambien 
porque carecen de poder predictivo en general. 

Una aproximacion alternativa es formular modelos mas simples, menos 
realistas pero manteniendo los rasgos fundamentals del sistema ffsico. 
Tales modelos estan descritos por unos pocos parametros, siendo por ello 
mas faciles de tratar y haciendo posible un estudio global del espacio de 
parametros sin perder intuicion fi'sica. 

Mediante la formulacion de estos modelos se pretende obtener resulta- 
dos relevantes sobre sistemas complejos, estudiando aproximaciones relati- 
vamente simples, como se suele hacer en Mecanica Estadi'stica (ME) donde 
sistemas tan complicados como redes ferromagneticas comparten muchas 
propiedades con modelos tan simples como el de Ising. En ME la relacion 
entre ambos sistemas se entiende a traves del Grupo de Renormalizacion. 
Desgraciadamente estamos lejos de poder probar una relacion de este tipo 
para modelos de flujo de partfculas, debido fundamentalmente a que no 
tienen un Ifmite termodinamico bien definido. Sin embargo la modelizacion 
de estos sistemas es necesaria. 

Los problemas que se encuentran cuando se estudian sistemas que en- 
vuelven flujo de partfculas estan relacionados con procesos de congestion. 
Los sistemas realizan una transicion de una situacion de trafico fluido a 
otra caracterizada por el atasco. 

Podrfamos utilizar las herramientas de la ME para estudiar estos prob- 
lemas si asimilamos este cambio a una transicion de fase. Un primer paso 
es buscar las variables que gobiernan la transicion. 

El flujo en la red esta condicionado por el numero maximo de partfculas 
que un nodo es capaz de contener, por la interaccion entre las partfculas y 
por la geometrfa de la red (por ejemplo el numero de coordinacion). 
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En nuestro modelo hay una inyeccion continua de particulas en la red, 
para cada una de ellas se elije al azar un nodo que sera su destino final. 
El movimiento de las particulas se implementa permitiendoles moverse en 
la red de acuerdo con deterniinadas reglas que definirenios mas tarde. En 
el trafico real las particulas encuentran a menudo obstaculos en su camino. 
Estos obstaculos se simulan limitando el niimero de particulas que un nodo 
de la red puede almacenar simultaneamente. 

Volviendo a la discusion desde el punto de vista de la ME, vamos a enfo- 
car la situacion en la que la densidad de particulas es alta. Intuitivamente 
es claro que una partfcula podrfa evitar zonas de atasco si le permitimos 
rodear el obstaculo. El sistema estara globalmente mas descongestionado y 
esperamos una mejora en el nrimero de particulas que alcanzan finalmente 
su destino. En terminos de ME, la discusion previa significa que calentando 
el sistema (T ^ 0) se obtendran mejores resultados 1^ 

En lo que concierne a la interaccion entre las particulas, henios intro- 
ducido un termino de contacto, la existencia de obstaculos en forma de 
nodos completamente ocupados. Las consecuencias de esta interaccion se 
pueden ver como si existiera un potencial infinito que actua en los nodos 
saturados. 

Un termino de interaccion apropiado se obtiene asignando carga electrica 
a cada partfcula. La repulsion alejara a las particulas de las regiones densa- 
mente cargadas. Consideraremos este tipo de repulsion actuando solo sobre 
primeros vecinos en la red. 

Este tipo de fuerzas repulsivas se consideran a menudo en problemas de 
polfmeros 1^. 

Hemos introducido un modelo en el cual hay un termino que controla 
la magnitud de las fluctuaciones termicas y otro que Simula una fuerza 
repulsiva. Los parametros que controlan estos terminos son respectivamente 
la temperatura y la carga. Veremos como ajustando ambos parametros es 
posible mejorar el flujo de particulas. 
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Abstract 



We study a particle flow model which may be used to get insight into 
various real traffic problems. The model is implemented using a discrete 
lattice, in which particles move towards their destination, fluctuating about 
the minimal distance path. A repulsive interaction between particles is 
introduced so as to avoid the appearance of a traffic jam. We have studied 
the parameter space finding regions of fluid traffic, and saturated ones, both 
regions being separated by abrupt changes. The improvement of the system 
performance is also explored by introducing a non-constant potential acting 
on the particles. Finally, we deal with the behaviour of the system when 
temporary failures in transmission occur. 

5.2 The Model 

We consider a two-dimensional lattice with coordination number 4, and 
periodic boundary conditions. The particles live in the lattice sites, la- 
beled by n = (jiQ , ni ) and can move from its site to one of its four nearest 
neighbors at every time step. The maximum number of particles that a 
site can contain will be denoted by B, and will be kept fixed for all the 
simulation at a value B = b. We denote the occupation number of the site 
n by a[n). In this notation, the particles are prevented from moving to 
sites with a{n) = B. In addition we consider an input queue at each site 
for the particles waiting to be injected on the lattice. 

By analogy with statistical mechanics systems, we work with the inverse 
of the temperature, (3. We denote the charge of the particles by k and the 
probability of particle injection to the lattice by p. 

The dynamic of the system is as follows: 

1. A lattice site is chosen at random. 

2. A particle is added, with fixed probability p, to the queue of the site, 
waiting to be introduced to the lattice [| 

3. If the queue of the site is not empty, and ^{n) < B, a new particle is 
introduced in the lattice, and an endpoint assigned to it at random. 

4. All particles at the considered site, try to move towards one of its 
four neighbors. For a given particle located at the position n we 

^The size of the queue is as big as necessary, so that it contains all particles waiting 
to be fed at that lattice site 
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must assign a probability of it jumping to each of its neighbouring 
locations. This probability is given by: 

P(±/i) = 7Vexp(±^sign(n^^ - 71^) - kct(7i^)) , (5.1) 

where rt^ here signifies the /i coordinate of site n, is the endpoint 
of the particle being considered, and N is the normalization constant. 
To choose between possible destinations we use a Heat Bath al- 
gorithm. 

The factor multiplying /3, is a potential term. It implies a constant 
force acting on the particle driving it to its endpoint. The n term pro- 
duces a repulsion between particles sitting in nearer neighbors sites. 
Obviously there is a wide range of potentials that could be considered, 
in order to produce more effective forces, and partial improvements 
will be expected. 

5. Movement is allowed if the chosen site has cr{n) < B, otherwise the 
particle remains at its original site until a new movement is attempted 
in the next iteration. 

6. A particle is removed from the lattice when it reaches its endpoint. 

The smaller p is, the weaker is the effect of the interaction between 
the particles. For fixed (/3, k) the fluid flow of particles on the lattice only 
takes place for those values of p being less than a certain threshold which is 
{(3, k) dependent. Above this threshold, the density becomes too high and 
the transmission process is prevented. We say the system is saturated. The 
saturation mechanism begins with the appearance of saturated domains, 
which grows in size making the movement of the particles more and more 
difficult. 

Our purpose is to quantify this threshold density, as well as to describe 
the flow properties along the parameter space (/?, K,p). 

From the point of view of Telecommunication Networks, the lattice sites 
can be thought of as nodes of such networks, and particles here can be in- 
terpreted as messages /packets. Within this approach, B can be identified 
with the available buffer space at the considered node. In this way, our 
model can be considered as a first step towards an abstract modeling of 
packet /message Telecommunication Networks. In such networks, packets 
are routed by each node according to routing tables dynamically maintained 
by the network so as to minimize some cost function along the trajectory of 
the packets In the present model, this cost function is simply the trav- 
eled distance, and, as the network state does not change during simulations, 
routing tables need not be dynamically calculated. A natural extension of 
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this model would be to allow variable link lengths between the nodes and 
to compute routing tables according to the path lengths. The role of the 
temperature parameter is to quantify how firmly the network will try to 
stick to its policy of minimizing the cost along a trajectory: as will be 
shown below, implementing the routing according to the cost minimization 
procedure in a rigid way (zero temperature), although this is the best naive 
choice from an individual user point of view, can be detrimental to the 
global behavior of the network and thus to its collective utility (the total 
throughput offered by the network). 

Congestion is an unavoidable phenomenon in uncontrolled packet net- 
works, since packets introduced at one node have no guarantee of finding 
the necessary resources (available buffer space) in the transit nodes. If 
uncontrolled, congestion result in packet losses which, in the case of data 
transmission, are detected by upper layer protocols (the Transport Control 
Protocol, TCP, for instance in the case of the Internet) and are corrected by 
retransmision request. Hence these packet losses will trigger more packets 
to be introduced into the network, thus further amplifying the congestion 
Two approaches are possible to avoid this vicious circle: reducing the source 
emission rate when a congestion is detected (this is the TCP approach) or 
trying to control the congestion inside the network by preventing pack- 
ets from accessing overloaded areas (this approach does not rely on any 
source behavior: in the extreme congestion case the source is simply denied 
any access to the network). This last option, congestion control inside the 
network, is presently a very active area of research in the networking com- 
munity To implement such a control, the overload information should 
be transmitted by a node to its neighbors at least. This is precisely what 
is modeled by the repulsive charge term in the present model: access to 
a loaded node is discouraged, and access to an overloaded node is simply 
forbidden. It should be noted that in the model as it stands, a node is 
instantaneously aware of the load state of its neighbors. A more complex 
model should take into account the latency introduced by the node-to-node 
transit time, leading to a description of the congestion control by some kind 
of retarded potential. It will be particularly important to take this latency 
aspect into account for modeling high speed networks |10 , but this is left 
for future work. 

The system we are considering here is also related to a lattice gas with 
fixed number of particles in a constant electric field E |jll[] . The electric cur- 
rent is related to the traffic speed. Also in this model a transition appears, 
between a disordered phase and another one were particles move collec- 
tively in the E direction. In our case the electric field is no longer constant, 
it depends dynamically on the particle environment: the surrounding elec- 
trical charge and the end point associated with the particle. As is usual for 
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traffic models, the number of particles is not constant; as a consequence, 
the transition mechanism is different. 

5.3 Numerical Simulation 

We have performed Monte Carlo (MC) simulations in order to study the 
parameter space. Here we present results obtained using a L x L lattice 
with L — 32 and periodic boundary conditions. The computations have 
been carried out on workstations. 

The starting configuration is obtained by generating a particle with 
probability p at each site of the lattice, therefore, the total number of sites 
occupied initially is about p x L^. The random number generator is based 
on one described in The time step is identified with a MC iteration, 
that is the update oil? lattice sites in the way described in section 2. 

The temporal evolution of the system exhibits a transient regime, char- 
acterized by the instability in the observables (see next section). After this, 
the system falls into an extremely long-living metastable state, where the 
flow properties do not change significantly with the temporal evolution. We 
say the system has reached the asymptotic regime. 

The time the system spends in the transient regime depends on the 
parameter space point k,p). In this regime, in the non-saturated region, 
the particle density is initially low and grows with the particle injection. 

This transit time also grows near the parameter space points where 
an abrupt change in the properties is observed (e.g. near the threshold 
density), reaching in this case up to 2 x 10''. 

For each value of the parameters, we have performed typically 8 x 10^ 
MC iterations. For /3 and k G (0, 0.4) the transient regime takes around 
5000 iterations, while above 0.5 for both parameters, this time falls to 
400 — 600 iterations. We have also performed the simulations starting from 
different configurations, allowing the system to evolve for up to 2 x 10^ iter- 
ations. We have computed the errors by calculating the dispersion between 
the results obtained from starting from these different configurations. 

5.4 Observables 

A correct description of the system is obtained from the measurement of the 
relevant observables. From their temporal evolution we are able to evaluate 
when the asymptotic regime is reached, or even whether this regime will 
or will not be saturated. From the averaged value of these magnitudes we 
obtain a quantitative description of the flow process. 
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Figure 5.1: Temporal evolution of M bellow the critical injection Pc- 



For a given configuration, we define the occupation M as: 
M = (l/y)f^aW. 

i=l 

Where V is L x L. 

The statistical average over the configurations (labeled by j) is: 



(5.2) 



N 



(M)= hm (l/iV)^M,- . 



(5.3) 



A'' being the number of averaged configurations. 

During the time interval tj — U = At, particles will reach their 
endpoint. We define the Band Width {Bw) as the number of particles that 
arrive at the endpoint per time unit. 

T) 

Bwm = ^ . (5.4) 

The statistical average for Bw is obtained from its mean value over a num- 
ber of time intervals Nt: 



{Bw)= lim (I/ATt) VBv^(j) . 

JVt— ><x) 



(5.5) 
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Figure 5.2: Temporal evolution of M above Pc- 



We can obtain a good measure of the system performance from the mean 
time taken by the particles to reache their endpoint, Tm- We compute Tm 
in the time interval At by adding the individual time spent by each particle 
(delay time), divided by n/: 



TM(At) 



(5.6) 



In the same way we define the statistical average as over Nt intervals as: 



{Tm) 



lim (l/iVr) VTm(j) 



(5.7) 



The occupation frequency of a certain occupation number, a{n), is de- 
fined as the number of times that the occupation a{n) appears at any lattice 
site. 



5.5 Phase Diagram 

We examine the parameter space (/3, K,p) searching for regions where sharp 
changes in the temporal evolution arise. At each (/?, k) value there is a p 
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Figure 5.3: Sections through the phase space cube. On the left side is 
plotted Pc versus P, on the right one, Pc versus k. 



value denoted by Pc such that for p < pc the asymptotic regime presents a 
stationary flow, and for p > pc the asymptotic regime is saturated and the 
flow is no longer possible. We plot the temporal evolution of M below pc 
(Figure 5.1 ) and above Pc (Figure 5.2) for some values of the parameters. 

This change is similar to a phase transition. The temporal evolution 
leads the system to one or another phase depending on the parameter space 
point. Once in the asymptotic regime, the non-saturated phase exhibits a 
dynamic equilibrium: the number of injected particles, equals the number 
of those arriving at their end position. This is reflected by (Bw) = px L^. 
In this phase, (Ta/) is constant, as well as (Af ) which is always less than B. 

The parameter space is divided into two regions by the surface defined 



0, 



by (/?, K,Pc)- Figure 5.3 shows two sections, for fixed (3 and k respectively. 

Above the surface, after the transient regime (M) = B and {Bw 
(Tm) diverges. 

It is possible to give a simple interpretation of this congestion phe- 
nomenon; when the load increases so does the probability that two particles 
residing in two neighboring saturated sites will want to exchange their sites, 
but of course they cannot do so since one of the particles would have to 
move first and then have nowhere to go to as no resource is will be available 
at the desired site. This mutual blocking then tends to propagate to other 
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neighboring sites and once above a given load a complete congestion will 
develop at the scale of the network. 

Such mutual blocking is well-known in most no loss routing schemes 
such as wormhole routing and is called deadlock ||l3|l . Deadlock avoidance 
is a very active research area, specially for massively parallel system inter- 
connections Q. 



P dependence 

The thermal fluctuations move particles away from their minimal distance 
path. The higher /3 is, the less the importance of these fluctuations. 



As shown in Figure 5.3 (left side), there is an interval of P values where 
the reduction of fluctuations have a positive influence on the throughput: 
Pc rises, as well as (Bw), while (Tm) decreases, as shown in Figure |3.4 

When going to higher /? values, the situation does not persist. The ab- 
sence of thermal fluctuations damages the throughput because particles are 
not able of going round obstacles. As a consequence Pc and {Bw) decrease. 
{Tm) does not appreciably change from its minimal value, corresponding 
to the infinite directionality one. 

Figure shows how the thermal fluctuations influence the saturation 
mechanism. If they are important, M grows slowly until it reaches B. If 
they are not significant a sharp jump appears in the temporal evolution of 
M, between its value in the transient regime and B. We conclude from this 
that the saturated domains grow faster in the absence of fluctuations, as 
would be expected from the earlier discussions. 

In particular, it is intuitive that a mutual deadlock will last longer if 
thermal fluctuations are disallowed (particles in mutual deadlock will re- 
peatedly attempt to exchange their sites) hence creating a larger local con- 
gestion area which can eventually evolve towards a global congestion. This 
was clearly observed in Q where only one buffer was available per site 
and this behavior is also exhibited in this model, hence illustrating the 
well-known fact [ p^ that deadlock is a consequence of our no loss routing 
scheme ^ and not a consequence of insufficient resources. 



K dependence 



The dependence of pc on k presents two different regions (see Figure 5^ ) . 

In the first one, Pc rises with increasing k. For these k values, the 
inclusion of the repulsion term helps the system to avoid congested regions. 
We denote by Kopt the value at which pc reaches its maximal value. 



movement is granted if and only if the available resource is available at the low end 
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Figure 5.5: Kopt) line. The errors are of the size of the k step measured, 
{Ak = 0.5.) 
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In the second region, k > Kopt, Pc decreases with k. Above Kopt on, the 
repulsion is too strong and the particles move far away from their minimal 
paths. As a result, the collapse appears for smaller injection densities. 

In Figure 5.5 we give the values of Kopt for some (3 values. The more a 
particle is restricted to its minimal path, the stronger the repulsion needed 
to avoid obstacles in the lattice. 

Therefore, there is a maximal value for the particle injection supported 
by each (3 value, which is reached for k in a neighborhood of Kopt- 



5.6 General Behavior and Optimization 

Let us focus on the behaviour of the other significant observables for the 
traffic flow, when (3 and k have been tuned to obtain a maximal Pc- 



Behavior of Bw 

As we have already pointed out, B\y in the asymptotic regime is propor- 
tional to p by a factor L^, as it may be checked in figure 6. Hence, at each 
(/3, k) value the maximum for (Bw) is reached when p = p^. 

Above Pc, (Bw) remains constant during the transient regime, and 
sharply falls to zero after it. Bw is a good measure to determine when 
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Figure 5.7: (Im) dependence on k. 
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Figure 5.8: Distribution of the occupation number a{n) times the frequency 
of this state versus a{n) after 5 x 10^ MC iterations. Transient regime 
contributions has been discarded. 
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Figure 5.9: /3 dependence of Pc for the distance dependent force (dashed 
line). Dotted line represents the values obtained for the constant force. 



the system has reached the asymptotic regime. 



Behavior of Tm 

The greater p the stronger the effect of the interaction over the particles: 
they move further and further away from their minimal path with increasing 
(Tm) (Figure 5_^). In Figure 5/7, we observe how (Tm) slightly increases 
with K until Hopt is reached From Kopt on, (Tm) increases with a high slope. 

We conclude that better performances are obtained for the maximal 
injection supported and for the bandwidth on the line {(3, Kopt)- The re- 
pulsion term damages Tj^i: the greater k, the larger the time the particles 
take to reach their destination. However, it is only for n > Hopt that (Tj^i) 
increases in a dramatic fashion. 



Figure shows the particle distribution in the simulation for /3 = 0.2 
and f3 — 4.0. The plot of a{n)heq[a{n)] exhibits a maximum around (M), 
and also reveals a wider distribution of particles for small /3 values. 
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Figure 5.10: < Tm > versus for fixed k = 1 in p w for the distance 
dependent force (dashed line) and constant force (dotted hne). 



5.7 Improving the throughput: distance de- 
pendent force. 

Although our implementation of the physical system has been inspired by 
the search for the simplest model accomplishing the desired features, par- 
tial improvements arc to be expected by taking into account details not 
included up to this point. As an example, we develop a possible improve- 
ment: quenching thermal fluctuations when they are no longer useful, that 
is, at the neighborhood of the endpoints. 

In the primary implementation, we have used a constant force and there- 
fore the particles support thermal fluctuations of the same strength no mat- 
ter what the remaining distance to the cndpoint is. We can expect a globally 
less congested lattice if particles which are only a few sites away from their 
destination are prevented from fluctuating. We will see how the inclusion 
of this feature does not spoil the good properties of the throughput, and 
also that the general behaviour of the system is not altered. 

A possible implementation to take this fact into account is obtained 
by introducing a distance dependence in the probability distribution, such 
that contributions of the fluctuations decrease with decreasing distance to 
the endpoint. 
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Figure 5.11: Delay distribution for the whole particles after 4 x 10^ MC 

iterations. Graphics on the left (constant force) must be compared with 
the right side ones (distance dependent force). The distributions in the top 
are calculated in (/3 = 0.2, k= l,p = 0.08) in the bottom they correspond 
to (/3 = 4,K= 1,}5 = 0.008). 
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A straightforward way of doing this is to include a dependence on the 
relative distance to the endpoint: the probability distribution now reads: 



P(±/i) = 7Vexp(±/3sign(7i^'^ - n^) ~ Ka{n^)) ( ^ ) ■ (5.8) 
Where r„ is defined by: 

'^(«o,ni) = \/ {n^o - "-o)^ + {nii- mY , (5.9) 



In Figure 5.9 the evolution of pc with (i is shown. A global throughput 
improvement is reflected by higher pc values. This effect is more remarkable 
when the size of the thermal fluctuations is important (small (3 values), 
corroborating our first intuition on the effect of fluctuations in the steps 
preceding the endpoint. 



In Figure |5.10| we plot the (3 dependence of Tm- We observe a global 
decrease in this time for all (3 values. 

The improvements conc erning the delay time, are not restricted to a 



smaller {Tm)- Figure 5.11 shows the delay distribution for all particles, 
compared with the delay distribution obtained with the constant potential. 
We see at (3 = 0.2 how the dispersion of the distribution strongly decreases 
when using a distance dependent force. So, the particles arrive in more 
similar times, increasing the uniformity and the reliability of the traffic 
flow. 



Shift 


n = 10 


n = 100 


n = 200 


13 = 4.0 K= 1.0 


ArPc 


0.09 


0.40 


0.58 


A,(Tm)(« Pc) 


-0.05 


-0.38 


-0.89 


Ar{M){^ Pc) 


0.00 


0.37 


0.55 


(3 = 0.2 n= 1.0 


ArPc 


0.04 


0.24 


0.41 


Ar {Tm){^ Pc) 


-0.00 


-0.08 


-0.37 


Ar{M){^ Pc) 


0.00 


0.43 


0.50 



Table 5.1: Shifted values of relevant observables for some values of the 
number of failures n. 



145 



0.02 0.04 0.06 0.08 0.100.02 0.04 0.06 0.08 0.10 




0.002 0.004 0.006 0.008 0.010 0.002 0.004 0.006 0.008 0.010 
P P 



Figure 5.12: (M) versus p for values of the fault percentage. Dotted lines 
represent (M) values obtained for an ideal lattice. On the top diagram the 
parameter space point is (/? = 4, k = 1), and on the bottom frames the 
plotted point is {(3 = 0.2, k ~ 1). 



5.8 Fault tolerance 

As we have pointed out, the lattice sites may simulate the control nodes of 
an information flow system. Needless to say the sites here are idealizations 
of nodes in real system because the possibility of communication failures is 
not allowed. Real nodes are subject to external factors, such as technical 
constraints or outside influences, that often damage the communication 
ability between some nodes. These nodes are then temporarily out of order, 
and cannot communicate or receive information from other nodes. It is 
therefore of interest to have an estimate of how robust a system is when 
temporary problems in the transmission occur. 

To implement the ocurrence of communication problems, information 
exchange is prevented at n randomly chosen sites during an interval of 
At = 10 MC iterations in the simulation. After this time interval these n 
sites are again allowed to communicate and another n sites are broken at 
random ^. 

In general, the measured values for observables will be shifted by an 



'in this section the probability distribution (5.1) is used 
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amount depending on n and in the parameter region. For each p value, we 
define the shift in the observable O as its value relative to the one obtained 
in the ideal system: 

We have studied the influence on the relevant parameters of communi- 
cation failures for n = 10, 100 and 200 with f3 = 0.2 and f3 = 4.0. 



Figure 5.12 shows the comparative evolution of (RI) for n = 10 and 
n = 100. 

The first consequence of a bad transmission is that the system must 
support globally higher occupation numbers. We measure smaller A,.pc 
values for lower /3, because, as has been already demonstrated, the lattice 
supports higher occupation numbers when the fluctuations are important. 

Ar{M) is almost zero for both /3 values with n = 10 (around 1% of 
nodes out of order), while with n — 100 and 200 Ar{M) is always large. 

In Table 5.1 we give the results obtained. As would be expected (Ta/) 



increases for all n values even though it is almost zero for n = 10. 



5.9 Summary and Outlook 

We have studied a model useful for describing the relevant processes occur- 
ring in traffic flow systems with immediate applications to network message 
passing and traffic problems in general. 

The introduction of the parameters /3 and k as controllers of the system 
behavior allows us to go a step further from purely descriptive models, 
because we are able to give prescriptions to improve the performance of the 
flow process. 

Deeper studies are also possible. Concretely an accurate study of the 
scaling with L of the relevant magnitudes as well as a detailed description 
of how the saturation time behaves, could lead us to the definition of quan- 
tities analogous to critical exponents. Also non-equilibrium states could 
be studied in order to monitor the parameters controlling the saturation 
process. 

The introduction of the k parameter has implied that the particles are 
able to avoid congested regions. The study has been limited to short-ranged 
interactions, because the particles at a site only see the occupation of its 
nearest neighbors. By informing the particles about the occupation of wider 
surrounding regions improvements in the throughput are expected. 
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Conclusiones 



Las conclusiones de esta memoria se han expuesto individualmente a lo 
largo de cada capitulo. A continuacion pasamos a resumir y comentar los 
resultados mas importantcs que se han obtenido. 

• El Modelo 0(4) Anti-ferromagnetico 

La introduccion de un acoplo negativo a segundos vecinos en el modelo 
a no lineal hace que aparezca frustracion en los vacios encontrados en 
el diagrama dc fascs. La consccucncia dc la frustracion cn cstc modelo 
es la aparicion de antiferromagnetismo en distintas dimensionalidades. 
La forma que hemos encontrado mas simple para introducir antiferro- 
magnetismo cn 4(i es trabajar en una red _F4 con acoplos negativos. La 
transicion del vacio desordenado paramagnetico al vacio AF a pianos 
es la mas interesante puesto que los datos numericos son compati- 
bles con segundo orden y los exponentes criticos no son los de campo 
medio. 

Se han expuesto las explicaciones alternativas, razonables dentro de la 
evidcncia numcrica, que pucdcii dar cucnta dc los resultados numericos 
Probablemente la explicacion menos plausible es la posibilidad de 
tener trivialidad logaritmica puesto que la dimension anomala del 
campo que hemos encontrado cs dcmasiado distinta dc ccro, y cstc 
tipo de correcciones dan cuenta de pequenas desviaciones a las predic- 
ciones de campo medio. 

Como sc ha senalado, a la vista de los resultados obtenidos en 0(2) 

y 0(3) cs muy dificil mantcner esperanzas de que la transicion en 
0(4) pueda ser de segundo orden. El escenario de primer orden debil 
parece ser el mas plausible. 

• La transicion de fase en el modelo SU(2)-Higgs 

Esta transicion dc fasc dc nucvo cs un ejemplo del problcma que se 
encuentra en d = 4 a la hora de discernir entre transiciones de primer 
orden muy debiles y transiciones continuas. 



Hemos visto que en SU(2)-Higgs la pnieba mas directa que dcmuestra 
el caracter de primer orden de la transicion, a saber la medicion del 
calor latente, no es accesible a un tiempo de calculo razonable. Se ha 
estudiado el problema en un espacio dc parametros extendido, lo cual 
ha ayudado a tener una vision global del mecanismo de debilitamiento 
de la transicion de fase. La observacion de la tendencia de los expo- 
nentes cn'ticos efectivos hacia los valores de primer orden constituye 
una evidencia del caracter de primer orden de la transicion. 

• U(l) compacto con topologias toroidal y esferica 

Este es un problema podriamos decir que clasico en el estudio de TCC 
en la red. Las cuestiones abiertas son dos: el orden de la transicion de 
fase y los mecanismos que la producen. En concreto se ha venido con- 
jeturando sobre una posible influencia de las condiciones de contorno 
en el orden de la transicion de fase. 

En esta memoria se ha expuesto el estudio de la transicion de fase 
desconfinantc cn U(l) puro gauge cn redes con topologia toroidal y 
esferica (homotopica a 5^). Se han estudiado los efectos de tamaho 
finito asociados a ambas topologias encontrando que en la red esferica 
hay mas efectos dc tamaiio finito incontrolados que en la toroidal 
dcbido a las inhomogcncidades. 

En cuanto al orden de la transicion de fase, se ha chequeado la exis- 
tencia de gap de energia para redes suficientemente grandes en ambas 
topologias, lo dial caracteriza a la transicion como de primer orden. 
Asi mismo se ha dado una estimacion del calor latente en el Ifmite de 
volumen infinito. 

• Modelos bidimensionales de flujo de particulas 

Se ha formulado un modelo para describir, en general, sistemas que 
envuelven un flujo de informacion. 

Se ha encontrado que se pueden obtener mejoras en el flujo permi- 
tiendo a las particulas fluctuar alrededor de la trayectoria de minima 
distancia. La introduccion del parametro k permite al sistema evitar 

regiones altamente congestionadas. 

La combinacion apropiada de flexibilidad en la direccionalidad del 
movimiento de las particulas hacia su destino, es decir fluctuaciones 

tcrmicas, y posibilidad dc evitar zonas congestionadas. cs decir el 
parrhetro k, es lo que permite un flujo optimo de particulas en la red. 
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evicogo 



En efecto, rematado ya su juicio, vino a dar en el mas 
extrano pensamiento que jamas dio loco en el mundo, 
y fue que le parecio convenible y necesario, asi para 
el aumento de su honra como para el servicio de su 
republica, hacerse caballero andante e irse por todo 
el mundo con sus armas y su caballo a buscar aven- 
turas y a ejercitarse en todo el aquello que el habia 
leido que los caballeros andantes se ejercitaban, des- 
faciendo todo genero de agravio, y poniendose en oca- 
siones y peligros, donde, acabandolos, cobrase eterno 
nombre y fama. 

Capitulo primero: Que trata de la condicion y cjcrcicio del famoso hidalgo 

Don Quijote de la Mancha 
"El Ingenioso hidalgo Don Quijote de la Mancha" 

(Miguel de Cervantes Saavedra) 
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